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AUTHOR’S  SUMMARY 

i 

The  aim  of  this  study  is  to  examine  the  strength  behaviour  of  an  important  con¬ 
structional  element,  the  lug,  when  fatigue  cracked,  and  to  propose  a  stress  intensity 
solution  considering  the  various  lug  geometries. 

First  the  stress  intensity  is  described  on  the  basis  of  the  linear-elastic 
principles,  which  is  followed  by  showing  their  significance  for  the  strength  behaviour  in 
fatigue  cracked  condition.  The  investigation  of  various  methods  for  determining  the 
stress  intensity  shows  the  crack  propagation  method  to  be  particularly  appropriate  for  the 
lug. 

Due  to  a  lack  of  suitable  tes  s  on  cracked  lugs,  a  gradual  approach  to  the  lug  is 
made,  starting  with  an  infinite  plate  with  a  hole  and  continuing  with  a  finite  plate  with 
and  without  pin  loaded  holes,  since  stress  intensity  solutions  are  already  available  in 
this  case. 

An  experimental  investigation  of  the  lug  not  only  enables  the  stress  intensity  to 
be  determined,  but  also  shows  the  interrelationship  between  the  essential  lug  parameters, 
the  crack  propagation  behaviour  and  the  stress  distribution  in  a  non-cracked  lug. 

Finally,  it  is  possible  to  evolve  a  stress  intensity  solution  for  cracked  lugs  in  the 
form  of  a  semi -empirical  approach,  which  takes  into  account  the  essential  lug  parameters. 
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GENERAL 

Fracture  mechanics,  ie  calculation  of  crack  propagation  and  the  critical  crack  length 
are  gaining  increasingly  in  importance  in  the  most  varied  technical  areas.  The  reasons  are 
to  be  found  primarily  in  stricter  safety  requirements  and  economic  considerations.  While, 
for  instance,  in  aircraft,  reactor  and  pipeline  construction  the  immediate  safety  of  man 
and  environment  is  of  primary  impoi Lance,  it  is  noticeable  that  fracture  mechanics  are 
applied  for  economic  reasons  in  areas  such  as  rolling  mills,  drilling  installations,  etc. 

A  shutdown  of  a  large-scale  plant  due  to  the  unforeseen  failure  of  a  component 
can  involve  a  company  in  considerable  losses  per  day. 

Basically,  fracture  mechanics  can  be  applied  wherever  damage  development,  ie  crack 
propagation  due  to  dynamic  loading  or  through  stress  crack  corrosion  under  static  load, 
cannot  be  excluded. 

Here  it  is  the  task  of  fracture  mechanics  to  predict  development  of  damage  and 
its  critical  quantity  under  service  conditions  in  order  to  be  able  to  take  suitable 
precautions . 

A  marked  drive  in  the  further  development  of  applied  fracture  mechanics  is  also 
likely  to  come  from  the  area  of  aircraft  construction*.  Fracture  mechanical  character¬ 
istics  have  been  used  for  a  long  time  as  a  criterion  in  the  selection  of  suitable 
materials. 

Inspection  intervals  for  aircraft  components  are  frequently  derived  from  crack 
2 

propagation  calculations  .  In  the  United  States,  Airforce  requirements  have  already  been 
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defined  in  Mil-Spec  83444  by  which  aircraft  manufacturers  must  prove  that  assuming  a 
crack  a  component  will  not  fail  under  service  conditions  within  a  certain  period.  Com¬ 
ponents  which  cannot  be  inspected  must  actually  tolerate  twofold  fatigue  life  with  crack. 

Due  to  these  requirements  fracture  mechanics  has  become  a  new  dimensioning  factor  in 
the  USAF  in  addition  to  strength  and  conventional  life  calculations,  the  latter  based  on 
Woehler  curves  and  linear  damage  accumulation  as  Miner.  Damage  tolerant  structure  can  be 
achieved  by  'fail  safe'  (redundancy)  or  'safe  crack  growth'  requirements.  Safe  crack 
growth  is  achieved  structurally  both  by  a  low  stress  level  and  crack  inhibitors  and  by 
suitable  selection  of  materials. 

The  following  paper  is  intended  to  provide  a  contribution  towards  the  improvement 
of  life  prediction  for  a  typical  aircraft  and  machine  element  within  the  framework  of 
damage  tolerant  design  philosophy. 

INTRODUCTION 

To  solve  structural  and  static  problems  the  lug  is  used  by  the  designer  in  many 
ways  as  an  element  for  transmitting  tensile  and  compressive  forces,  bu’t  without  bending. 

4-8 

There  are  various  methods  of  calculating  the  fatigue  strength  of  this  structural 
element.  The  known  literature  does  not  at  present  include  a  stress  intensity  solution 
for  different  lug  geometries,  assuming  a  unilateral  crack  with  straight  crack  front. 
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In  practice,  the  cracks  occurring  in  lugs  are  predominantly  quarter-elliptic  edge 
cracks  or  semi-elliptic  surface  cracks.  If  in  the  light  of  dimensi 'tiing  philosophy  (damage 
tolerance)  it  is  accepted  that  cracks  may  be  present  from  the  start,  then  the  probability 
of  the  presence  of  one  crack  in  the  critical  lug  cross-section  is  greater  than  that  of  a 
crack  on  both  sides.  For  this  reason  it  is  the  aim  of  this  paper  to  develop  a  method  of 
calculation  for  the  unilaterally  cracked  lug. 

Starting  with  the  stress  intensity  solution  for  a  unilateral  crack  with  straight 
linear  crack  front,  more  complicated  crack  forms  could  be  considered  like  the  edge  or 
surface  crack  for  instance  by  linear  superposition  with  known  answers  for  cracks  with 
curved  crack  fronts.  Regarding  a  crack  with  straight  or  at  least  approximately  straight 
crack  front  has  the  advantange  of  reducing  the  investigation  to  an  essentially  two- 
dimensional  problem. 

1  FUNDAMENTALS  OF  LINEAR  ELASTIC  FRACTURE  MECHANICS 

9  1 0 

Apart  from  the  material  characteristics  ’  which  will  not  be  gone  into  here,  the 
stress  intensity  plays  a  decisive  part  in  the  calculation  of  crack  propagation  and 
residual  strength  or  critical  crack  length.  It  forms  the  basis  of  every  fracture  mechanics 
calculation.  The  theoretical  fundamentals  of  the  stress  intensity  within  the  frame¬ 

work  of  linear  elastic  fracture  mechanics  are  outlined  below. 

The  whole  derivation  of  the  elastic  stress  field  at  the  crack  tip  has  already 
been  compiled  in  Ref  11.  The  following  sunmary  has  been  extracted  from  Ref  12. 

1 . 1  The  Airy  stress  function 

Within  a  coordinate  system,  X,  Y,  Z  of  a  stressed  body  the  stresses 


xy 


T 

yz 


can  be  defined  at  any  point  (x,y,z). 


For  the  plane  stress  state 


a  =  t  =  t  =  0  . 

z  xz  yz 


For  the  plane  strain  state  e  =  0  ,  a 

z  z 


v(a  +  o  )  . 
x  y 


For  the  plane  stress  state  the  following  equilibrium  conditions  result  (see  Fig  1.2) 
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If  the  displacements  in  x  and  y  direction  are  designated  u  and  v  ,  the  result  for 
the  strains  is  (see  Fig  1.3), 
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The  equations  between  stress  and  strain  in  the  plane  stress  state  are  formulated 
according  to  Hooke's  Law  by 
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Equation  (1-1)  is  satisfied  when 
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Function  C  is  termed  Airy's  stress  function. 

If  equations  (1-2)  and  (1-4)  are  inserted  in  (1-3)  and  differentiated  twice  the 
result  is  the  compatibility  equation: 
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(1-6) 


In  general  a  plane  linear  elastic  problem  can  be  solved  by  finding  the  stress  function 
which  satisfies  equation  (1-6).  In  addition  the  stresses  calculated  from  (1-4)  must 
satisfy  the  boundary  conditions  of  the  problem.  The  stress  function  for  a  special  case 
must  be  estimated  on  the  basis  of  experience.  Approximation  solutions  are  described  in 
detail  in  Ref  13. 

1 . 2  Complex  stress  function 

A  complex  stress  function  can  be  defined  as  follows: 


Z(z)  =  ReZ+ilmZ;  z  =  x  +  iy  .  (1-7) 

Z  should  be  an  analytical  function  with  a  unique  derivative  dZ/dz  .  This  leads  to  the 
Cauchy-Riemann  conditions: 
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3ReZ  =  3ImZ  dZ 

3x  3y  6  dz 
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(1-8) 


Various  complex  forms  of  the  Airy  stress  function  can  be  used  to  solve  crack  problems 

In  the  case  of  crack  type  I  (see  Fig  1.1)  the  function  suggested  by  Westergaard' 
18  19 

is  suitable.  While  Sih  and  Eftis  and  Liebowitz  demonstrated  that  the  Westergaard 
function  is  not  exact,  this  has  no  effect  as  far  as  the  singular  expressions  of  the 
stress  are  concerned. 

The  Westergaard  function  reads: 


i|>  =  ReZ  +  ylmZ 


(1-9) 


Z,  Z  and  Z'  being  defined  by: 


d!  -  dZ_  dZ 

dz  ’  dz  ~  ’  dz 


0-10) 


With  the  Cauchy-Riemann  equations  (1-8)  it  follows 

V2ReZ  =  V2ImZ  =  0  . 


(1-11) 


This  means  that  (1-9)  automatically  fulfils  the  compatibility  equation  (1-6). 
Using  (1-4)  produces  for  stresses 
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(1-12) 


Every  analytical  function  Z(z)  will  lead  to  stresses  as  in  equations  (1-12).  What  is 
required,  however,  is  to  find  a  function  Z(z)  for  the  problem  in  question  which  also 
fulfils  the  boundary  conditions. 

1 . 3  Solving  the  crack  problems 

For  an  infinite  plate  under  biaxial  stress  the  following  stress  function  is 
obtained  for  crack  type  1  : 
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,  z  .  . 

Z  =  —  -  ,  z  =  x  +  ly 

>J -  a^ 


(1-13) 


The  function  is  apart  from  (-  a  <  x  <  a,  y  =  0)  analytical.  The  edge  stresses  follow 

from  (1-12).  At  infinity  where  I Z I  -+  °°  ,  o  =  o  =  a  and  x  =  0  .  At  the  crack 

x  y  xy 

face  a  =  t  =  0  .  This  means  that  the  boundary  conditions  are  satisfied, 
y  xy 

In  order  not  to  have  to  replace  Z  by  (z  +  a)  it  is  logical  to  shift  the  co¬ 
ordinate  system  to  the  crack  tip.  In  the  general  case,  ie  without  taking  the  boundary 
conditions  into  account,  Z  must  take  the  following  form: 


Z  = 


f  (z) 
/z 


(1-14) 


f(z)  must  be  real  and  constant  at  the  origin. 

In  accordance  with  (1-12)  a  =  t  =0  applies  at  the  crack  surface. 

y  xy 

The  required  real  term  and  the  constant  value  of  f(z)  is  given  at  the  crack  tip 
by  Kj  ,  since 


z  bo 


/2ttz 


0-13) 


i  0 

Using  polar  coordinates  at  the  origin  with  Z  =  Re  ,  the  stresses  at  the  crack  tip  can 
be  calculated  from  equations  ( 1  - 1 2)  and  (1-15). 
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Parameter  in  the  above  equations  is  termed  stress  intensity  factor.  For  r  -*  0  in 

the  elastic  case  the  stresses  become  infinite.  The  stress  intensity  factor  is  therefore 
a  measure  of  the  stress  singularity  at  the  crack  tip.  Since  the  material  behaviour  is 
assumed  to  be  linear  elastic  the  stress  nust  be  proportional  to  the  external  loading.  In 
the  case  of  uniaxial  load  with  j  in  the  infinite  Kj  must  be  proportional  to  a  . 

Kj  must  also  be  proportional  to  the  root  of  the  crack  length.  For  an  infinite  plate  the 
sole  characteristic  length  is  the  size  of  the  crack.  Consequently  must  have  the 

following  form: 

K[  =  c-»a  .  (1-17) 
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In  the  case  of  biaxial  load  the  stress  function  is  given  by  (1—13) .  The  shift  in  the 
origin  of  the  coordinate  system  to  the  crack  tip  changes  (1-13)  to 


»  z  (  z  +  2a  ) 


Comparison  of  (1-15)  with  (1-18)  shows  that 


K  =  oi’a 


(l-!r>) 


As  one  may  assume  that  the  stress  system  parallel  to  the  crack  is  not  disturbed  by  the 
crack,  the  solution  for  uniaxial  loading  must  be  the  same  as  for  biaxial  loading.  In 
the  case  of  uniaxial  loading  factor  C  in  (1-17)  assumes  the  value  »r"  .  In  addition  to 
the  stresses,  the  displacements  for  the  plane  strain  state  can  also  be  determined  from 
equations  (1-2)  and  (1-12). 


This  leads  to 
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Equations  (1-16)  are  exact  solutions  for  the  stress  field  in  the  range  r  551  0  .  These 
solutions  can  still  be  used  in  a  range  in  which  r  is  small  compared  with  the  crack 
length. 

2  STRESS  INTENSITY  AS  BASIS  FOR  CALCULATION  OF  RESIDUAL  STRENGTH  AND  CRACK 
PROGRESSION 

2. 1  Unstable  crack  propagation 

Calculation  of  the  fracture  behaviour  of  cracked  components  was  carried  out  bv 
Griffith22. 

The  calculation  rests  on  the  assumption  that  fracture  occurs  in  ideally  brittle 
materials  when  the  elastic  energy  released  by  crack  extension  is  exactly  the  same  as  the 
energy  needed  for  crack  propagation.  This  energy  formula  can  be  described  as  follows. 
Imagine  an  infinite  plate  with  a  crack  going  through  the  centre.  The  plate  is  loaded 
with  a  uniaxial  stress  a  applied  at  infinity.  The  potential  energy  U  of  the  entire 
system  can  then  be  described  by  the  formula 
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U  =  Un  -  U  +  U 
0  a  y 


where  =  elastic  energy  of  the  non-cracked  plate 


U  =  reduction  of  the  elastic  energy  bv  the  insertion  of  a  crack 
a 

U  =  increase  in  surface  energy  due  to  enlargement  of  crack. 

.23  . 

With  the  aid  of  a  stress  calculation  developed  by  Inglis  ,  Griffith  was  able 


to  show  that 


Moreover  the  surface  energy  is  equal  to  the  product  of  the  specific  surface  energy 

of  the  material  y  and  the  new  surface  of  the  crack 
e 

U  =  2(2ay  )  .  (2-3 

Y  e 

The  entire  elastic  energy  of  the  system  follows  from  this. 

U  =  U0-^L-4aYe  •  ^ 

Derivation  of  3U/3a  produces  the  equilibrium  conditions  for  crack  extension 


Crack  extension  in  ideally  brittle  materials  is  therefore  determined  by  the  product  of 

the  external  stress,  the  root  of  the  crack  length  and  by  material  characteristics. 

With  E  and  y  as  material  characteristics  the  right-hand  side  of  the  equation  (2-5) 
e 

corresponds  to  a  constant  which  is  characteristic  of  a  specified  ideally  brittle  material. 
It  follows  from  this  that  crack  extension  in  these  materials  occurs  when  the  product  of 
a/a  reaches  a  critical  value.  This  value  can  be  determined  experimentally  for  instance 
with  a  plate  with  central  crack  by  measuring  a  and  a  at  the  moment  of  crack  extension 
and  transferred  to  other  specimen  geometries,  taking  account  of  the  effect  of  the  geometry. 
It  was  only  the  latter  which  made  practical  application  of  fracture  mechanics  possible  in 
a  wide  variety  of  technical  areas. 

Equation  (2-5)  can  be  rearranged  as  follows: 


the  left-hand  side  being  termed  energy  release  rate  G  .  The  right-hand  side  of  the 
equation  represents  the  crack  resistance  of  the  material. 

In  1948  Irwin"*4  suggested  that  Griffith's  fracture  criterion  for  ideally  brittle 
materials  should  be  transferred  by  modification  to  brittle  materials  and  to  metals 
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exhibitin'"  plastic 
Orowan  at  the  same 
is  composed  of  the 
sequently  equation 


deformations  at  the  crack  tip.  A  similar  suggestion  was  made  by 
time.  According  to  this  the  crack  resistance  of  technical  materials 
sum  of  the  surface  energy  and  the  plastic  deformation  action.  Con- 
(2-6)  was  changed  accordingly. 


VE 


»'2(y„ 


(2-7) 


A  crack  becomes  unstable  when  the  left-hand  side  of  the  equation  (2-7)  reaches  a  critical 
value,  ic  the  critical  stress  intensity  K^,  .  For  crack  type  I  the  critical  stress 
intensity  is  designated  below  with  .  For  the  infinite  plate  with  central  crack  2  a 

and  external  stress  3  fracture  occurs  when  equation 


cj/Ta  =  K 


1C 


(2-8) 


is  satisfied.  Ref  25  describes  the  experimental  method  of  determining  Che  critical  stress 
intensity  .  It  is  usual  in  the  literature  to  designate  the  critical  stress  intensity 
in  plane  strain  state  K ^  and  in  plane  stress  state  . 

Crack  type  I  which  occurs  most  frequently  in  practice  is  always  used  as  basis. 


The  critical  stress  intensity  in  the  range  of  the  plane  strain  state  is 

primarily  a  material  constant.  In  the  range  of  the  plane  stress  state,  however,  the 
critical  stress  intensity  depends  on  the  material  thickness. 


For  the  critical  stress  intensity  in  plane  stress  and  in  plane  strain  state  it  is 
assumed  that  the  plastic  zone  in  front  of  the  crack  tip  is  small  compared  with  the 
length  of  the  crack  and  the  dimensions  of  the  component.  This  is  shown  in  the  test 
by  the  linear  relationship  between  force  and  crack  expansion.  If  linearity  no  longer 
exists  (e lasto-plastic)  then  the  critical  stress  intensity  is  dependent  on  additional 
factors.  Relevant  tests  were  carried  out  in  Ref  26. 
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Critical  stress  intensity  values  are  compiled  in  Ref  27  for  a  variety  of  materials. 

As  soon  as  one  changes  from  the  infinite  plate  with  central  crack  to  specimens 
resembling  components  the  effect  of  the  crack  form,  crack  arrangement,  specimen  geometry 
and  type  of  external  load  must  be  accounted  for  by  a  correction  function  Y  .  Correction 
functions  for  the  various  crack  models  are  compiled  inRefs  28  and  29. 

K  =  o/iraY  .  (2-9) 


2 . 2  Stable  crack  propagation 

Paris  ’  demonstrated  the  connection  between  stable  crack  propagation  due  to 
cyclical  loading  and  stress  intensity. 

By  means  of  crack  propagation  tests  with  constant  stress  amplitude  Paris 
determined  the  crack  growth  per  load  cycle.  In  double  logarithmic  paper  plotting  the 
pr.oagation  rate  against  the  related  stress  intensity  factor  range  in  the  middle 
section  B  of  the  curve  shown  in  Fig  2.2  produced  a  linear  relation.  Deviations  only 
occur  with  very  small  AK^  values  in  section  A  and  in  the  vicinity  of  the  instability 
point  in  section  C.  This  enabled  the  technically  interesting  area  to  be  described  by 
a  simple  exponential  function: 


da 

dN 


C  (AK  )np 
P  I 


(2-10) 


It  has  been  shown  for  a  number  of  high  strength  alloys  that  the  effect  of  the  stress 
amplitude  in  the  range  of  fluctuating  tensile  stresses  on  crack  propagation  behaviour 
is  adequately  reflected  by  the  Paris  equation  for  a  defined  stress  ratio  R  . 

The  material  characteristics  C  and  n  in  the  Paris  equation  are  determined  by 

P  P 

means  of  constant  amplitude  crack  propagation  tests.  These  characteristics  depend 
essentially  on  the  material  and  stress  ratio  R  . 

32  .  . 

Forman  expanded  this  equation  m  that  he  took  account  of  the  effect  of  the  mean 
stress  and  the  crack  propagation  behaviour  on  reaching  the  critical  crack  length: 


da 

dN 


np 

VAV 

O  -  R)KIC  -  AKX 


(2-11) 


The  Forman  equation  is  used  very  frequently  in  practice.  MBB  set  up  a  programme  to 

33 

determine  the  Forman  characteristics  .  This  programme  has  been  included  in  the  structure 
calculation  handbook  (HSB)  of  the  German  aviation  industry. 

In  practice,  however,  loading  sequences  occur  frequently  with  variable  amplitudes 

34-39 

and  mean  stress.  It  has  been  demonstrated  m  numerous  tests  that  crack  propagation 

on  transition  from  high  to  lower  stress  peaks  can  be  delayed  appreciably.  For  this 
.  40  4 1 

reason,  eg  Willenborg  and  Wheeler  have  developed  crack  propagation  retardation  models 


14 

for  the  variable  load  sequence.  Besides  file  Forman  equation  these  retardation  models  have 

been  evaluated  by  means  of  extensive  test  results  which  are  typical  for  aircraft  con- 
4!’ 

struction  and  suggestions  worked  out  for  calculating  the  crack  propagation  under 
operational  loads.  The  retardation  models  are  also  based  on  the  stress  intensity. 

3  METHODS  OF  DETERMINING  STRESS  INTENSITY 

As  already  shown  in  sections  2.1  and  2.2,  knowledge  of  the  variation  of  stress 
intensity  is  of  considerable  importance  to  the  calculation  of  crack  propagation  and 
residual  strength.  Extensive  investigations  have  therefore  been  carried  out  to  determine 
the  stress  intensity  for  various  components  or  specimens  resembling  components.  Various 
methods  were  used.  Analytical  methods  can  be  used  for  very  simple  geometries.  However, 
in  the  case  of  geometries  with  complex  boundary  conditions  numerical  and  experimental 
methods  are  frequently  employed. 

Various  methods  of  determining  the  stress  intensity  are  presented  below  and 
evaluated  in  regard  to  their  applicability  to  the  problem  'lug'.  Additional  information 
on  several  methods  is  contained  in  Ref  43. 

3. 1  Analytical  and  numerical  methods 

3.1.1  Analytical  methods 

The  elementary  equations  to  describe  the  stress  and  displacement  field  at  the 
crack  tip  for  the  case  "infinite  plate  with  central  crack"  wore  sot  up  witli  the  aid  of 
these  methods. 

These  results  still  serve  as  the  starting  point  for  the  solution  of  various  crack 
prob lems . 

The  fact  that  the  stress  and  displacement  field  for  any  crack  problem  always 
assumes  the  same  shape  offers  the  opportunity  of  determining  the  stress  intensity 
i  ndi rec  1 1 y  . 

From  the  point  of  view  of  the  engineer,  however,  the  analytical  methods  are  of 
secondary  importance.  in  general  the  houndary  conditions  are  detected  precisely  by 
analytical  methods.  However,  this  is  usually  only  possible  for  an  infinite  plate  with 
central  crack..  In  determining  stress  intensity  by  means  ol  analytical  methods  the 
attempt  is  made  to  tied  the  Airy  stress  function  corresponding  to  the  crack  problem. 

For  crack  type  I  the  stress  function  of  Westergaard *  is  frequently  used, 

=  ReZ  '  vlmZ  (3-1) 

The  solution  of  this  function  was  given  in  section  1.  Muskhel ishvi 1 i ' " ,  among  others, 
suggests  a  complex  stress  function 


’1  =  Rp|Zij>(z)  +  J'Mz)dzl 


(3-2) 


Using  (1-4)  it  follows 


a  +  o  -  4  Re  If'  (  z )  |  . 

x  y 


(  3-1) 
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Thus  it  can  be  shown  that 


2K. 


lim  (o  +  ci  ) 
n  x  v 
r-*-0 


0 

COS  77 
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This  method  was  tested  in  detail  by  Sih  and  used  by  Erdogan 


(3- A) 


If  a  simple  solution  of  these  equations  is  not  possible  numerical  methods 
(difference  equations)  can  be  used  as  approximation  solutions. 

46 

To  solve  the  technically  interesting  problem  "crack  from  hole"  Bowie  used  the 
method  of  conformal  mapping. 

3.1.2  Method  of  finite  elements  as  numerical  process 

The  method  of  finite  elements  has  been  used  to  an  increasing  extent  recently  to 

...  47 

calculate  the  stress  field  in  the  vicinity  of  the  crack  tip  .  The  versatility  of  this 
method  enables  complicated  component  geometries  to  be  calculated.  The  fundamentals  and 
applications  of  the  method  of  finite  elements  are  described  in  Ref  48. 


In  the  method  of  finite  elements  the  elastic  continuum  is  replaced  by  a  limited 
number  of  elements  which  are  only  connected  at  their  nodal  points. 

Forces  or  displacements  between  the  elements  can  only  be  transmitted  through 
these  nodal  points. 

The  displacements  of  the  nodes  are  the  unknown  quantities  in  an  appropriate  equation 
system. 


In  the  plane  stress  state  the  sole  displacements  are  u  and  v  .  In  the  simplest 
case  it  is  assumed  that  these  displacements  are  linear,  eg  u  =  ax  +  by  +  c  , 
v  =  ex  +  fy  +  g  . 


In  many  cases  triangular  elements  are  used  with  a  node  at  each  angle  of  the  element. 
The  displacements  of  the  nodal  points  are  Uj,  U^,  and  Vf  ,  V^,  . 

The  displacements  of  the  nodal  points  must  satisfy  the  equations  for  u  and  v  . 

By  inserting  the  node  displacements  and  node  coordinates  six  equations  are  obtained  with 
six  unknowns  a  to  g  for  an  element. 

Solving  the  equation  system  produces  all  the  displacements  within  the  element. 

These  assumptions  ensure  continuity  between  neighbouring  elements,  ie  the  linear  change 
in  displacements  along  the  contact  line  of  two  elements  will  be  the  same  in  both  elements 
because  of  the  common  nodal  points.  It  follows  from  ex  =  3U/3x  and  e y  =  3l’/3y  that 
with  linear  displacements  deformations  are  constant  within  an  element.  Of  course 
functions  of  a  higher  order  are  possible  to  describe  the  deformations  within  each 
element.  This  considerably  improves  the  accuracy  of  the  FE  calculation,  especially  in 
the  area  of  steep  stress  gradients. 


The  bonding  forces  between  the  elements  can  also  be  presented  as  nodal  displacements. 
Finally,  equilibrium  equations  have  to  be  set  up  for  the  nodal  points. 


16 


In  a  plane  problem  the  nodal  forces  have  two  components,  namely  in  x  and  y 
direction.  The  equilibrium  equations  are  yielded  by  the  sum  of  the  component  forces 
of  various  elements  which  meet  at  the  node. 

The  forces  occurring  at  the  nodes  of  the  marginal  elements  are  equated  with  the 
external  loads  or  stresses. 

The  resulting  equation  system  can  be  solved  by  electronic  computers. 

Basically  there  are  two  different  approximations  to  determine  the  stress  intensity 
with  the  aid  of  the  method  of  finite  elements. 

There  is  the  direct  method  in  which  the  stress  intensity  is  calculated  from  the 
stress  or  displacement  field. 

In  the  indirect  method  the  stress  intensity  is  determined  via  the  relationship 
with  other  parameters  such  as  stiffness,  elastic  energy  or  J  integral. 

For  the  direct  method  the  result  of  general  analytical  solutions  of  crack  problems 
is  used.  The  stress  intensity  can  thus  be  calculated  from  the  stresses  and  displacements 
in  the  vicinity  of  the  crack  tip  with  the  aid  of  equations  (1-19)  and  (1-21): 


=  a 


ij  fTTe) 


u. 

1 

c/rf.  (0) 


(3-5) 


The  necessary  stresses  and/or  displacements  are  obtained  from  calculation  by  the 
finite  element  method.  The  stress  for  an  element  in  the  vicinity  of  the  crack  tip,  the 
distance  r  of  the  centre  of  gravity  of  the  element  from  the  crack  tip  and  the  angle  6 
(see  Fig  1.2)  are  inserted  in  equation  (3-5),  thus  yielding  the  stress  intensity  at  a 
certain  crack  length  a  . 

The  stress  intensity  can  be  calculated  in  a  similar  manner  from  the  displacement. 

In  the  ideal  case  the  stresses  and  displacements  of  different  elements  should  lead 
to  the  same  Kj.  value,  provided  the  elements  are  in  the  area  of  validity  of  the  analyti¬ 
cal  solution  (3-5),  ie  in  the  immediate  vicinity  of  the  crack  tip. 

To  reduce  the  integration  error  owing  to  the  finite  dimensions  of  the  finite 
elements  very  fine  subdivision  of  the  grid  in  the  crack  area  is  recommended  if  elements 
with  linear  displacement  are  used. 

To  model  the  crack  tip  singularity  special  crack  tip  elements  were  developed  from 

the  so-called  isoparametric  elements.  For  instance  the  FE  programme  developed  by 
49 

MARC  contains  a  rectangular  element  with  eight  nodes  and  nine  integration  points 
within  the  element. 

Singularity  is  approximated  by  shifting  together  two  angle  nodes  situated  on 
one  side. 

The  advantage  of  these  elements  is  to  be  seen  in  that  far  fewer  elements  are  needed 
for  idealisation  of  the  crack  tip  environment  with  the  same  accuracy  than  is  the  case  for 
elements  with  linear  di  placement  formula. 

I _ _ 
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For  this,  approximately  6-8  of  these  elements  are  usually  required  at  the  crack  tip. 
The  computing  costs  can  also  be  reduced  substantially  by  this  means. 

7he  stress  intensity  can  also  be  calculated  indirectly  via  changes  in  stiffness 
with  increasing  crack  length.  The  procedure  corresponds  to  the  experimental  compliance 
method  described  below.  Here  the  test  to  be  performed  in  sect 'on  3.2.1  is  replaced  by 
the  FE  calculation. 

The  FE  calculation  provides  the  displacement  of  the  force  introduction  points  for 
various  crack  lengths. 

The  change  in  stiffness  as  a  function  of  crack  length  is  thus  obtained,  as  in 
Fig  3.7. 

The  K  value  is  calculated  to  equation  (3-6). 

The  J  integral  can  also  be  used  for  indirect  calculation  of  the  K.  value. 

2  /  1 

As  long  as  elastic  deformations  are  concerned  the  equation  /E  =  Cj  =  J  applies. 

For  example  Chan,  et  al~>®  determined  the  stress  intensity  with  the  J  integral 
for  a  compact  tensile  specimen.  The  external  contour  of  the  specimen  was  used  as 
integration  path. 

Strain  energy  density  was  determined  from  the  nodal  point  stresses  and  the  nodal 
point  forces  used  as  surface  line  forces. 

The  indirect  methods  have  the  advantage  that  extrapolation  is  not  required  as  in 
the  direct  method  and  according  to  the  known  literature  no  particularly  dense  element 
grid  is  said  to  be  necessary. 

3. 2  Experimental  methods  of  determining  stress  intensity 

3.2.1  The  compliance  method 

In  the  compliance  method  the  stiffness  of  a  component  at  various  crack  lengths  is 
measured  experimentally.  The  stress  intensity  is  calculated  via  the  crack  extension 
force.  According  to  Irwin  the  elastic  energy  U  changes  by  the  value  on  extension 

of  the  crack.  This  value  is  termed  crack  extension  force  by  Irwin: 

311 


From  the  force  P  and  the  displacement  AS  of  two  points,  generally  the  force 
introduction  points,  the  elastic  energy  U  is  derived. 


PAs 

o 


The  displacement  As  can  be  replaced  by  the  compliance 


r" 1  1  _  As 

C  P 


(3-7) 
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Equations  (3-7)  and  (3-8)  produce  for 


u 


e 


If  (3-9)  is  inserted  in  (3-6)  it  follows 


2 

P  dC 


da 


=  C. 


(3-9) 


(3-10) 


Via  the  known  relationship  between  stress  intensity  K  and  crack  extension  force  G 


GI  -  t 


plane  stress  state 


(3-10)  can  be  converted  to 
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0^  =  —  (I  -  v  )  plane  strain  state 
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The  schematic  diagrams  A,  B  and  C  in  Figs  3  to  7  show  the  procedure  in  experi¬ 
mental  determination  of  the  stress  intensity  . 

Fig  A  shows  the  specimen  extension  As  with  increasing  force  p  at  different 
crack  lengths  a  .  With  increasing  length  of  crack  the  stiffness  of  the  specimen 
decreases,  ie  the  specimen  or  component  becomes  'softer'. 

In  B  the  reciprocal  values  of  the  slopes  from  A  are  plotted  over  the  length 
of  crack  a  .  A  distinct  increase  in  compliance  C  *  with  increasing  crack  length  a 
can  be  discerned  in  this  diagram  too.  If  one  differentiates  in  B  compliance  C  1 
with  respect  to  crack  length  a  ,  ze  forms  dC  '/da  and  plots  this  in  C  over  the 
relevant  crack  length,  then  the  stress  intensity  K  can  be  calculated  with  this  curve 
and  with  equation  (3-11). 

In  the  main  the  stress  intensity  can  be  determined  only  for  cracks  with  linear 
crack  front  with  the  compliance  method. 

For  semi-elliptical  surface  or  quarter-elliptical  edge  cracks  the  compliance 
method  would  only  yield  a  median  value  for  the  entire  crack  front. 
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3.2.2  Determination  of  stress  intensity  with  crack  propagation  measurements 

It  has  been  demonstrated  in  numerous  constant  amplitude  crack  propagation  tests 
that  the  crack  propagation  rate  can  be  presented  as  a  function  of  the  stress  intensity 


This  means  that  the  same  crack  propagation  rate  is  always  measured,  irrespective  of 
crack  length,  stress  and  specimen  geometry,  for  a  specific  stress  intensity  range.  If 
the  stress  intensity  is  to  be  determined  for  a  component,  a  crack  propagation  test 
under  constant  amplitude  loading  must  first  be  performed  with  the  component.  From  the 
crack  propagation  curve  the  crack  propagation  rate  is  determined  as  a  function  of  the 
crack  length. 

A  constant  amplitude  crack  propagation  test  is  carried  out  with  a  test  specimen 

which  agrees  with  the  component  in  respect  of  material  and  thickness  under  the  same 

conditions  (stress  ratio,  loading  frequency,  environment).  Since  the  stress  intensity 
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for  the  specimen  must  be  known  the  crack  propagation  rate  can  be  presented  from  the 
crack  propagation  curve  as  a  function  of  the  stress  intensity  range, 


da/dN  =  f(AK  ) 


(3-12) 


This  crack  propagation  characteristic,  however,  can  be  transferred  to  the  component, 
as  stated  at  the  beginning. 

If  now  a  certain  value  for  the  crack  propagation  rate  has  been  determined  in  the 
component  at  a  certain  crack  length,  then  this  must  have  been  brought  about  by  a  specific 
stress  intensity  range  according  to  (3-12).  Therefore  the  appropriate  stress  intensity 
range  can  be  allotted  to  the  crack  propagation  rate  established  on  the  component. 

Fig  3.8  shows  the  procedure  in  the  form  of  a  flow  diagram. 

It  is  customary  in  the  literature  to  quote  the  correction  function  as  non- 
dimensional  quantity  depending  on  crack  length  instead  of  the  stress  intensity.  The 
stress  intensity  is  used  as  reference  quantity  which  would  result  at  the  same  stress 
and  crack  length  in  an  infinite  plate  with  central  crack  under  uniaxial  lead. 


AKIMKI0 
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(  3-14) 


As  in  the  Compliance  method,  this  process  is  suitable  primarily  for  cracks  with  linear 
crack  front. 

This  has  an  important  advantage  in  that  this  process  provides  crack  propagation 
curves  for  the  component  investigated  with  which  the  accuracy  of  the  stress  intensity 
or  correction  function  determined  can  be  verified  by  re-checking  the  figures. 
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3.2.3  Photoelasticity  and  strain  gauges 

With  the  aid  of  equations  (1-19)  and  (1-21)  the  stress  intensity  can  be  calculated 
by  knowing  the  stresses  or  strains  in  the  area  of  the  crack  tip.  Stresses  or  strains  can 
be  determined  experimentally  by  means  of  photoelasticity  or  with  the  aid  of  strain 
gauges. 

In  contrast  to  the  above  experimental  methods,  photoelasticity  offers  the  possibil- 

5  3— 5  6 

ity  of  dealing  with  cracks  with  curved  fronts  by  'stress  freezing' 

In  this  method  the  crack  is  replaced  by  a  notch  with  a  finite  radius.  The  differ¬ 
ence  from  the  fatigue  crack  is  covered  by  appropriate  correction  factors. 

The  stress  intensity  is  determined  with  the  aid  of  the  shear  stresses.  The  shear 
stresses  themselves  are  obtained  from  the  isochromatics  (=  lines  of  constant  principal 
stress  differences)  and  the  related  fringe  number. 


K 


I 


(3-15) 


Since  the  value  cannot  be  determined  in  the  immediate  vicinity  of  the  notch 

(r  =  0),  the  K  values  must  first  be  calculated  at  some  distance  from  the  notch.  A  good 
approximation  value  for  K  at  r  =  0  is  obtained  by  extrapolation  (see  Fig  3.10). 


! 

i 


Fundamentally  any  method  of  determining  stress  intensity  can  be  used  which  enables 
stresses  or  strains  to  be  measured  near  the  crack  tip. 
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For  the  strain  gauge  method  several  strain  gauges  are  glued  in  front  of  the 

crack  tip  and  strains  e  and  e  measured.  The  stresses  are  yielded  by  the  following 

x  y 

equations : 

a  =  — ^  (e  +  ve  )  ,  a  =  — tt  (c  +  vc  )  .  (3-16) 

x  2  x  v  v  ,2yx 

l-i  J  J  1  -  v  1 


With  these  stresses  the  stress  intensity  can  be  determined  in  the  same  way  as  by  the 
photoelastic  method. 

With  this  method  care  must  be  taken  to  ensure  that  the  strain  gauges  lie  outside 
the  plastic  zone  of  the  crack  tip. 

4  INVESTIGATIONS  KNOWN  FROM  THE  LITERATURE  FOR  DETERMINING  STRESS  INTENSITY  FOR 
CRACKS  EMANATING  FROM  THE  EDGE  OF  THE  HOLE 


In  the  most  varied  technical  areas  fatigue  cracks  in  dynamically  stressed  structures 
have  frequently  proved  to  emanate  from  holes.  According  to  an  investigation  of  aircraft 
frames  of  the  USAF^,  30%  of  all  fatigue  cracks  discovered  occurred  in  holes  loaded  with 
pins  or  rivets.  Causes  of  crack  generation  are  usually  the  relatively  high  stress 
concentration,  the  generation  of  fretting  corrosion  in  loaded  holes  and  faulty  production 
of  the  hole.  Appropriate  stress  intensity  solutions  were  therefore  determined  relatively 
early  for  holes  affected  by  cracks,  compared  with  other  configurations  of  cracks  of 
interest  in  practice. 
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This  applies  to  the  infinite  plate  as  well  as  the  finite  plate  with  hole,  solutions 
being  available  both  for  the  loaded  and  unloaded  hole. 

The  lug  can  be  regarded  as  a  special  case  of  a  cracked  hole.  This  is  due  to  the 
finite  dimensions  in  width  and  head  height  and  the  introduction  of  force  which  occurs 
exclusively  by  means  of  a  pin. 

Since  lug  junctions  in  structures  are  frequently  not  redundant,  the  failure  of  this 
force  transmitting  element  can  cause  failure  of  the  entire  structure.  In  aircraft, 
failure  of  a  lug  can  even  lead  to  a  crash. 

Although  the  lug  in  dynamically  loaded  structures  to  a  certain  extent  represents 
an  increased  risk  and  is  thus  of  particular  interest  for  the  application  of  fracture 
mechanics,  there  is  so  far  no  appropriate  stress  intensity  solution  in  the  known 
literature.  Therefore  it  appears  sensible  to  approach  the  more  complex  problem  of  the 
lug  by  starting  from  simpler  similar  crack  configurations  with  known  solutions.  The 
effect  of  the  hole,  finite  dimensions  and  loaded  hole  on  the  stress  intensity  will  be 
illustrated  step  by  step. 

4. 1  Infinite  plate  with  unloaded  hole 

For  the  case  of  an  unloaded  hole  with  uni-  and  bilateral  crack  in  an  infinite 
plate  Bowie^*  determined  the  stress  intensity  analytically  by  conformal  mapping  and 
using  the  Airy  stress  function: 

Kj  =  i'vVaV  .  (4-1) 

The  correction  function  Y  contains  the  crack  length  ratio  a/R  as  parameter.  The 
crack  emanating  from  the  edge  of  the  hole  is  designated  a  and  the  radius  of  the  hole  R 
(see  Fig  4.2). 

Y  =  f (a/R)  .  (4-2) 

The  correction  function  for  the  unilateral  crack  Y^  and  the  bilateral  crack  Y,  was 
indicated  by  Bowie  in  the  form  of  tabular  values  (see  Table  4.1). 

The  Bowie  correction  can  be  expressed  by  the  following  equations  with  an  accuracy 
of  i2Z  (see  Ref  63) : 


Y j  =  0.707  -  0.187  +  6.55V2  -  10.54V3  +  6.85V4 


Y,,  =  1.0  -  0.15V  +  3 . 46  v  “  -  4.47V  +  3.52V 


(4-3) 

(4-4) 


V 


(4-5) 


If  Table  4.1  is  studied  it  will  be  seen  that  close  to  the  edge  of  the  hole  the 
correction  function  compared  with  that  of  a  central  crack  of  the  same  length  is  also 
greater  in  an  infinite  plate. 
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With  a  crack  length  ratio  of 


a/r  -  2 


the  influence  of  the  hole  on  the  stress  intensity  can  be  ignored  to  a  very  great  extent. 

4. 2  Infinite  plate  with  hole  and  introduction  of  force  by  means  of  a  pin 

•  64 

In  most  cases  holes  are  loaded  with  rivet  or  screw  connections.  Grandt 
developed  a  stress  intensity  solution  for  the  'loaded  hole'  and  compared  this  with  the 
Bowie  solution  for  open  holes. 


In  determining  the  stress  intensity  Grandt  used  an  analytical  process  suggested 
by  Rice^"’  which  is  briefly  described  below. 

4.2.1  Linear  superposition 

On  the  assumption  that  the  stress  intensity  Kj.  and  the  crack  opening  displacement 
are  known  for  a  particular  crack  geometry  and  particular  introduction  of  force,  the  stress 
intensity  can  be  calculated  for  any  required  introduction  of  force,  as  demonstrated  by 


According  to  Rice  the  stress  intensity  for  a  given  crack  geometry  under  any  load 
can  be  determined  to  the  following  formula: 


Kt  =  J'shdr  +  /  f  hdA  (4-6) 

IF  A 

s  is  the  stress  vector  which  at  the  edge  T  here  acts  on  the  crack  tip.  f  designates 
a  force  which  occurs  within  the  area  bordered  by  F  (see  Fig  4.3).  Vector  h  is  a 
universal  weight  function  for  the  given  crack  geometry  and  a  load  with  known  stress 
intensity.  The  weight  function  is  given  by 


For  the  constant  H 


h 


h(x,y ,a) 


H  3u 
2K];  3a  ’ 


(4-7) 


H 


E/(l 


for  the  plane  stress  state 

2 

v*")  for  the  plane  strain  state. 


The  modulus  of  elasticity  is  designated  E  and  the  Poisson  ratio  v  .  and  u  are 

the  known  stress  intensity  and  the  displacement  field  for  a  particular  load  configuration. 

4.2.2  Establishing  stress  intensity  with  the  aid  of  the  weight  function  for  the 
loaded  hole 

For  the  case  of  a  loaded  hole  equation  (4-6)  can  be  simplified  as  follows.  First 
the  stress  vector  s  is  determined.  In  doing  so  the  reaction  stress  p(x)  in  the 
future  crack  plane  is  used  instead  of  the  force  of  the  pin  acting  on  the  hole  radius. 

The  stress  distribution  p(x)  in  the  crack-free  state  can  be  calculated,  for  example,  by 
the  method  of  finite  elements. 
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The  stress  vector  obtained  is  then 


S  =  0 

x 


S  =<S  =  p  (x) 

i  y 


along  the  crack  planes 


S  =  S  =0  along  the  hole  circumference, 
x  y 


.(4-8) 


Since  no  volume  force  occurs  within  the  edge  F  ,  the  pin  force  being  replaced  by  the 
reaction  stress  p(x)  , 

f  =  0  .  (4-9) 


Since  the  displacement  field  u  shows  only  one  component  in  y  direction  due  to 
uniaxial  loading 

U  =  n  (4-10) 

y 

is  set  down. 

For  the  stress  intensity  of  the  loaded  hole 


a 


(4-11) 


is  obtained. 

Kg  is  the  stress  intensity  solution  of  Bowie  for  a  unilateral  crack  in  an 
unloaded  hole. 

n  is  the  appropriate  crack  opening  displacement.  Since  Bowie  did  not  give  data 
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on  the  crack  opening  displacement,  Grandt  used  an  approximat ion  : 


R  i 

n  -  4  -f  U/2n)! 


(4-12) 


t  is  defined  as 


t;  =  a  -  x  . 


(4-1  3) 


The  approximation  solution  as  equation  (4-12)  is  only  valid  within  the  range 

C  <  a/10  .  (4-14) 

The  stress  distribution  p(x)  was  calculated*^  for  the  loaded  hole  by  the  finite 
element  method. 

In  Fig  4.4  the  stress  intensity  is  shown  as  a  fun-tion  of  the  creek  length  for  the 
loaded  and  unloaded  hole  for  two  different  fits.  It  will  be  seen  that  with  increasing 
clearance  between  pin  and  hole  the  stress  intensity  also  increases. 


,1 


Fig  4.4  also  shows  that  up  to  about  a  crack  length  ratio  of  a/R  "  1  the  stress 
intensity  of  the  loaded  hole  is  greater  than  that  of  the  unloaded  hole.  With  increasing 
crack,  length,  on  the  other  hand,  the  stress  intensity  of  the  loaded  hole  drops  again 
below  that  of  the  unloaded  hole. 

In  contrast  to  previous  observations,  the  crack  propagation  direction  was  assumed 
under  an  angle  of  81°  based  on  the  FE  calculation,  and  not  perpendicular  to  the  loading 
direction  (see  Fig  4.4). 

4 . 3  Plates  with  finite  dimensions 

The  stress  intensity  solutions  of  Bowie  for  the  unloaded  hole  and  Grandt  for  the 
'loaded  hole'  can  only  be  used  if  the  hole  diameter  and  the  crack  are  small  compared 
with  the  dimensions  of  the  plate. 

If  this  is  not  the  case,  the  plate  edge  must  be  allowed  for  in  the  stress  intensity. 

This  applies  to  plates  with  loaded  and  unloaded  holes. 

4.3.1  Finite  plate  with  unloaded  hole 

Newman^  developed  a  stress  intensity  solution  for  a  plate  with  finite  width  and 
open  hole,  using  the  so-called  boundary  value  problem  method. 

The  method  is  based  on  finding  a  stress  function  which  the  biharmonic  equation  (.1-6) 
and  the  boundary  conditions  fulfil  in  a  number  of  points  along  the  edge  of  the  component 
(see  section  4. 3.2.1).  Thus  the  following  expression  of  the  stress  intensity  is 
obtained  for  the  uni-  or  bilateral  crack  in  an  unloaded  hole: 

ki  =  --Vb.Ac  •  (4_I5) 

The  correction  function  for  the  stress  intensity  allows  with 


Y 


W 


I  sec 


-  3R  +  a 
3  V  -  a 


(4-lb) 


for  the  finite  width  W  ,  with 


Y 


bl 


see  equation  (4-3) 


the  effect  ol  the  hole  according  to  Bowie  and  with  >  sec  (3R)/W  the  r<atio  between 
hole  diameter  3R  and  width  W  . 

Index  1  in  the  correction  function  designates  a  unilateral  and  Index  2  a  bilateral 

c  rack  . 

It  the  correction  function  of  Bowie 


Y 


B 


bl 


(see  equation  (4-3) 


is  compared  with  that  of  Newman 
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V  V  V  /  2Rtt 

YN  =  VBlV/seC  IT 

it  can  be  demonstrated  with  equations  (4-15)  and  (4-3)  that:  Y^  >  Y^  applies  over 
the  entire  crack  length  range. 

The  correction  function  and  thus  the  stress  intensity  according  to  Newman  becomes, 
compared  with  the  Bowie  solution,  the  greater,  the  greater  the  diameter  2R  in  relation 
to  width  W  and  the  nearer  the  crack  is  to  the  edge  of  the  plate. 

Fig  4.6  shows  for  2R  =  20  and  W  =  100  a  comparison  of  the  two  correction 
functions. 

4.3.2  Finite  plate  with  introduction  of  force  by  means  of  a  pin 

The  best  known  stress  intensity  solutions  for  finite  plates  with  introduction  of 

6  3 

force  by  a  pin  were  determined  by  Newman  by  the  boundary  value  problem  method  and 
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by  Broek  by  the  superposition  method  .  The  two  methods  will  be  described  briefly  below. 
4.3.2. 1  Boundary  value  problem  method*^’70 

This  is  a  numerical  method  which  is  used  to  determine  the  unknown  coefficients  in 
a  stress  function  series.  In  Ref  68  an  approximation  developed  by  Muskhelishvili  based 
on  complex  variables  was  used  for  the  stress  function  series.  With  the  aid  of  the  stress 
function  series  the  biharmonic  function 


=  0 


(4-17) 


in  which  equilibrium  and  compatibility  conditions  (see  section  1.1)  are  combined  can  be 
presented  analytically.  Due  to  symmetry  conditions  uneven  terms  of  the  stress  function 
series  can  be  eliminated.  The  progression  is  limited  to  a  certain  number  of  terras. 

The  coefficients  are  determined  through  the  appropriate  boundary  conditions.  The  boundary 
conditions  can  be  laid  down  in  the  form  of  stresses,  forces  or  displacements. 
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4. 3. 2. 2  Principle  of  superposition 

The  basis  of  the  principle  of  linear  superposition  is  that  the  stress  intensity 
solution  for  a  certain  crack  model  with  a  somewhat  more  complicated  form  of  loading  is 
built  up  from  already  known  solutions  for  the  same  crack  model  but  with  a  simpler  form 
of  loading. 

Fig  4.7  shows  this  process  with  a  finite  plate  with  hole  and  introduction  of 
force  by  means  of  a  pin. 

Thus  one  obtains  for  K. 

A 


Since 


ka  =  kb  +  kd-ke 


K  + 

B  D 


(4-18) 


2 


(4-19) 
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For  K_  the  Bowie  solution 

D 

known  solution 


(see  equations  (4-3)  and  (4-4)  can  be  used  and  for 


>  ’a 


the 


(4-20) 


The  boundary  effect  is  allowed  for  by  a  width  correction  Y(a  ,,/ W)  .  Thus  the  correction 

eft 

function  for  the  loaded  hole  has  the  following  form 


The  effective  crack  length  a  ^  's  defined  in  Fig  4.8. 

4.3. 3  Comparison  of  correction  function  for  the  finite  plate  with  hole,  with  and 
without  introduction  of  force 

The  correction  function  Y  of  Newman  and  of  Broek  was  determined  for  the  uni¬ 
lateral  crack  from  a  hole  in  a  finite  plate  with  force  introduction  via  a  pin.  Fig  4.9 
shows  that  Broek's  correction  function  yields  higher  values  for  Y  over  almost  the 
entire  crack  length  range  than  Newman's  solution.  The  same  tendency  was  seen  in  the 
unloaded  hole  (see  Fig  4.6). 

For  purposes  of  comparison  the  correction  function  was  given  for  the  bilateral 
crack  according  to  Cartwright^1  for  a  loaded  hole.  Cartwright  determined  the  correction 
function  by  the  Compliance  method  described  in  section  3.2.1.  For  short  crack  lengths 
Cartwright's  correction  function  agrees  well  with  that  of  Broek.  Only  at  greater  crack 
lengths  the  loss  of  supporting  effect  in  the  case  of  a  bilateral  crack  appears  to  make 
itself  felt  by  a  smaller  drop  in  the  correction  function. 

Fig  4.9  also  shows  the  correction  function  for  a  unilateral  crack  from  an 
unloaded  hole  in  a  finite  plate. 

Up  to  a  crack  length  ratio  of  '0.5  this  correction  function  is  distinctly  below 
that  of  a  loaded  hole.  A  reversal  takes  place  only  when  the  cracks  are  longer.  This 
was  also  observed  in  the  infinite  plate  with  loaded  and  unloaded  hole  (see  Fig  4.4). 

Irrespective  of  the  plate  edge  the  following  can  therefore  be  stated  as  generally 
applicable:  the  correction  function  Y  or  stress  intensity  for  the  loaded  hole 

is  greater  than  for  the  unloaded  hole  in  the  area  of  short  crack  lengths.  A  reversal 
takes  place  for  longer  crack  lengths. 

4 . 4  Lug 

Investigations  on  the  lug  which  have  only  recently  appeared  in  the  literature 
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have  been  undertaken  by  L.F.  Impellizerri  and  D.L.  Rich  “  and  A.F.  Liu  and  H.P.  Kan 

4.4.1  Investigation  by  Impellizor  i  and  Rich 

The  aim  of  the  investigation  was  to  carry  out  flight-by-flight  crack  propagation 
tests  on  a  lug  with  unilateral  crack  and  then  subject  the  test  result  to  a  mathematical 

check . 
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It  was  necessary  for  this  purpose  to  determine  the  stress  intensity  specifically 
for  the  lug  used  in  the  crack  propagation  test.  The  method  described  in  section  4.2.1 
was  used  for  the  purpose. 

The  function  derived  by  Bueckner  for  the  lateral  crack  was  used  as  the  weight 
function  (see  section  3.2.1). 

The  stress  distribution  was  determined  by  the  method  of  finite  elements  (see 
Fig  4.10). 

Data  on  the  idealisation  of  the  introduction  of  force  in  the  FE  calculation  are 
unfortunately  not  given.  As  shown  later,  however,  the  pattern  of  stress  and  thus  also 
the  stress  intensity  depend  on  just  this.  Fig  4.11  shows  the  pattern  of  the  stress 
intensity  or  correction  function  and  the  dimensions  of  the  lug. 

The  correction  function  includes  the  correction  factor  of  1.13  derived  by 

Bueckner  for  edge  cracks.  The  correction  factor  is  intended  to  allow  for  the  load-free 

surface  from  which  the  crack  emanates.  A  similar  factor  has  already  been  quoted  for 

'>8 

edge  cracks  by  Brown  and  Srawley  at  (1-12). 

The  correction  function  for  lugs  Y,  would  therefore  assume  the  value 

A 

Y,  -*  1.13  for  a  ->■  0  if  the  stress  intensity  were  based  on  the  stress  o  occurring 
A  •  max 

at  the  edge  of  the  hole. 

4.4.2  Investigation  of  Liu  and  Kan 

On  the  principle  of  linear  superposition  (see  section  4. 3. 2. 2)  Liu  and  Kan  first 
determined  the  correction  function  Y  ^  for  uni-  and  bilateral  cracks,  varying  internal 
diameter  in  relation  to  lug  width. 


Next  the  correction  function  was  in  addition  calculated  for  a  lug  with 
R/W  =  0.125  by  the  finite  element  method. 

Fig  4.13  shows  the  idealisation  of  the  lug  for  the  finite  element  calculation. 

With  the  aid  of  the  stress  intensity  determined  by  the  superposition  method  crack 
propagation  was  calculated  for  various  lugs  under  constant  amplitude  load  sequence  and 
compared  with  the  test  results. 

When  the  results  presented  in  Fig  4.12  are  regarded  two  facts  are  noted.  For  one 
thing  the  value  for  the  correction  function  Y^  for  a  -*■  0  is  considerably  higher 
than  that  determined  by  Impellizzeri  and  Rich. 


a  ->  0 


=  5  according  to  Impellizzeri  and  Rich 

8  -*■  1  4  according  to  Liu  and  Kan. 


The  investigations  by  Cartwright^'  and  Crandt^  already  mentioned  make  the  value 
determined  by  Impellizzeri  and  Rich  appear  more  realistic. 


The  reason  for  the  values  determined  by  Liu  and  Kan  for  Y  "  8  v  1 4  is  probably 

the  too  simple  approximation  process,  namely  linear  superposition.  The  complexity  of 

the  lug  seems  to  be  only  partially  covered  by  the  method  of  linear  superposition 
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already  used  by  Broek 


28 


It  is  probably  also  due  to  the  approximation  method  that  there  is  no  difference 
worth  mentioning  between  uni-  and  bilateral  cracks  in  regard  to  .  This  is  contra¬ 

dicted  by  results  by  Bowie  (see  Table  4.1)  as  also  Broek  and  Cartwright  (see  Fig  4.9). 
It  is,  however,  presupposed  that  the  difference  between  uni-  and  bilateral  cracks  in 
regard  to  stress  intensity  is  mainly  transferable  to  the  lug  from  infinite  plates  and 
strips.  In  view  of  the  loss  of  supporting  effect  consequent  on  bi-lateral  cracking, 
a  distinct  difference  ought  to  occur  with  increasing  crack  length  between  uni-  and 
bilateral  cracks  particularly  in  the  lug. 

The  results  of  Liu  and  Kan  show  furthermore  that  for  the  unilateral  as  for  the 
bilateral  crack  the  correction  function  Y 

A 


Y  -*•  <»  for 
A 


a 


W  -  2R 


2 


applies. 
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Crack  propagation  tests  on  lugs  have  shown  at  suitably  low  stress  that  in  a 

unilateral  crack,  crack  propagation  runs  steadily  through  one  cheek  of  the  lug  whereas 

the  opposite  undamaged  cheek  takes  over  more  and  more  load  until  it  finally  bears  the 

full  load.  This  leads  to  the  conclusion  that  for  a  lateral  crack  the  correction 

function  Y,  in  the  border  line  case 
A 

W  -  2R 


reaches  a  finite  value. 

What  is  surprising  about  the  results  of  Liu  and  Kan  is  that  for  a  lug  with 
R/W  =  0.125  checking  by  the  method  of  finite  elements  agrees  fairly  well  with  the  result 
of  the  superposition  method.  This  would  appear  to  be  purely  coincidental  since  the 
methods  used  here  are  fundamentally  different  approximation  solutions. 

If  one  assumes  that  by  using  a  crack  element  the  stress  intensity  can  be  achieved 

with  1-2Z  accuracy,  then  greater  deviations  can  be  attributed  largely  to  idealisation 
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of  the  introduction  of  force 

To  determine  the  stress  intensity  Liu  and  Kan  assumed  introduction  of  force  in 
the  form  of 

p  =  2P  cos(0/nrt)  .  (4-22) 

In  addition  a  concentrated  point  loading  introduction  of  force  as  well  as  one  of  the 
form  of 

p  =  3P  cos^(9/4rt)  (4-21) 


was  investigated. 


These  three  distributions  are  shown  in  Fig  4.14. 
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How  precisely  the  stress  intensity  was  determined  is  shown  finally  in  a  comparison 
carried  out  by  Liu  and  Kan  between  crack  propagation  test  and  calculation. 

The  crack  propagation  characteristics  used  for  the  mathematical  check  were  deter¬ 
mined  with  compact  tensile  specimens  (see  Fig  8.  ’). 

To  allow  for  the  material  scatter  occurring,  crack  propagation  for  the  lugs  was 
calculated  with  the  crack  propagation  characteristics  for  the  upper  and  lower  limit  of 
the  scatter  band. 

As  Fig  4.15  shows,  the  results  of  crack  propagation  tests  on  the  lug  lie  between 
the  upper  and  lower  limit  of  the  crack  propagation  curves  determined  mathematically. 

The  difference  between  the  e rack  propagation  calculation  with  crack  propagation 
characteristics  on  the  safe  and  unsafe  side  can  deviate  by  1-5  fold  from  the  test 
result.  These  are  erratic  values  which  completely  contradict  experience  hitherto  with 
constant  amplitude  crack  propagation  tests.  Deviations  ot  the  order  of  approximate  1 v 
10-20%,  on  the  other  hand,  would  be  realistic.  This  assumes,  howevei ,  that  the  specimens 
used  for  determination  of  the  crack  propagation  characteristics  are  taken  from  one  batch 
and  the  tests  performed  under  the  same  conditions. 

5  SUMMARY  OF  INVKSTIOATION'S  DRAWN  FROM  I'HK  L 1  TKRATl'KH 

A  table  of  the  crack  models  examined  and  methods  used  to  determine  the  stress 
intensity  is  contained  in  Fig  5.1. 

If  the  most  important  results  of  the  investigations  carried  out  in  Chapter  4  are 
summarised  purely  qualitatively,  the  following  el  lefts  can  be  identified  in  regard  to 
stress  intensity  and  correction  function  for  the  different  forms  of  plates  with  holes. 

The  ratio  between  length  of  crack  and  radius  of  hole  is  of  decisive  importance  tor 
the  correction  function.  The  greater  the  hole  radius,  the  greater  is  the  crack  length 
range  in  which  due  to  the  hole  the  correction  function  is  increased  compared  with  a  crack 
of  the  same  length  in  a  plate.  In  the  immediate  '’icinity  of  the  edge  of  the  hole,  tot- 
instance,  the  correction  function  for  an  unloaded  hole  with  unilateral  crack  in  an 
infinite  plate  is  for 

a  =  0  ;  Y(a/R)  =  1.14  (.see  Table  4 .  1  1  . 

With  increasing  distance  of  the  crack  tip  1 rom  the  edge  of  the  hole  the  effect  of  the 
hole  on  the  correction  function  also  decreases.  If  the  length  of  the  crack  amounts  to 
once  to  twice  the  hole  radius  the  effect  ot  the  hole  is  removed  to  a  great  extent. 

If  the  hole  is  loaded  with  a  pin  this  produces  an  additional  increase  in  the  stress 

intensity  or  correction  function  for  a  crack  in  the  immediate  vicinity  ot  the  edge  ot 
the  hole.  With  increasing  crack  length,  however,  the  correction  turn-tic'  is  again 
progressively  reduced  (see  Fig  4.4). 

'4  , 

According  to  Crandt  the  fit  be  tween  pin  and  hole  is  also  said  to  have  an  effect 

on  the  stress  intensity.  According  to  this  greater  clearance  between  pin  and  hole  would 

also  produce  higher  stress  intensity  (see  Fig  4.4).  This  mav  certainly  be  ot  importance 
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for  relatively  small  cracks.  For  longer  cracks,  however,  it  may  be  assumed  that  the  type 
of  fit  is  of  secondary  importance.  The  overestimation  of  the  effect  of  fit  by  Crandt  is 
chiefly  attributable  to  the  inadequacy  of  the  method  used  by  him  to  determine  the  stress 
intensity.  In  the  method  described  in  section  4.2.2  the  essential  quantity  used  is  the 
stress  distribution  in  a  crack-free  state  which  largely  depends  on  the  fit. 

A  further  effect  on  stress  intensity  arises  from  the  finite  size  of  the  plate, 
ie  from  the  edge  of  the  plate.  As  the  crack  tip  approaches  the  edge  of  the  plate 
the  stress  intensity  increases  again  (see  Fig  4.6). 

72  7  3 

According  to  the  investigations  available  in  the  literature  the  hole  diameter 

and  lug  width  are  to  be  regarded  as  essential  influencing  parameters  in  the  lug.  The 
head  height  of  the  lug  has  not  been  taken  into  account  as  a  parameter  in  any  investigation. 

Aside  from  the  preceding  purely  qualitative  description  of  the  stress  intensity  or 
correction  function  with  different  influencing  factors,  it  seems  possible,  i rrespec t i ve 
of  the  shape  and  loading  of  the  plate  with  hole  (see  Fig  5.1),  to  state  in  a  simple  manner 
quantitative  values  for  the  correction  function  at  the  edge  of  the  hole. 

If  the  investigation  of  Bowie  for  the  infinite  plate  with  hole  (see  section  4.1) 
and  of  Impellizzeri  for  the  lug  (see  section  4.4.1)  are  considered  again  the  following 
values  are  yielded  for  the  correction  functions: 


Plate  with 

hole: 

Y(a  = 

0) 

=  3.39 

according  to 

„  .  61 

Bowie 

Lug: 

Y(a  = 

0) 

=  5.65 

according  to 

Impel  1 izzeri ' ~ 

Furthermore  if  the  stress 

at  the 

edge 

of 

the  hole 

is  considered 

the  result  at  a 

stress  of  a,T  =  1  is  for 
N 

Plate  with 

hole : 

a 

- 

3.0 

according  to 

„  7s 

Petersen  and 

max 

3.06 

according  to 

N  e  ub  e  r  7 6 

Lug: 

a 

max 

= 

5.0 

according  to 

Impellizzeri' 

For  the  lug  the  stress  is  related  to  the  gross  cross-sort  ion. 


The  ratio  between  correction  function  and  stress  yields  in  both  cases,  . tor  the 
plate  with  hole  and  for  the  lug,  a  factor  of  1.1). 

For  plates  with  edge  or  surface  crack  this  laetor  is  intended  to  take  into  account 
that  the  edge  from  which  the  crack  emanates  is  load-free,  4.  no  force  occurs  parallel 
to  the  crack. 

On  the  other  hand,  in  an  infinite  plate  with  central  crack  stresses  occur 
parallel  to  the  crack  planes  which  suppress  the  opening  of  the  crack  to  a  minor  extent. 
The  stress  intensity  for  the  infinite  plate  with  central  crack  is  laid  down  by 


K 


1 


( 5-  1  ) 


If  the  centrally  cracked  plate  is  halved  by  sectioning  perpendicular  to  the  crack  plane 
the  result  is  two  plates  with  lateral  cracks  while  the  cut  edges  are  load-free.  The 
stress  intensity  here  is 

Kj  =  Yn/rTa  .  (5-3) 
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For  Y  values  are  quoted  between  1.12  and  1.13.  Neuber  for  instance  determined  a 
value  of  1.1215  for  Y  for  the  half  plate  with  crack  at  the  free  edge.  Since  the  Airy 
stress  function  was  determined  here  by  progressive  development  Y  should  be  regarded  as 
an  approximation  value  which  depends  to  a  certain  extent  on  the  number  of  terms  used. 
Bueckner^  quotes  this  value  as  1.13. 

If  the  nominal  stress  is  used,  which  can  be  related  to  either  the  gross  or 

net  cross-section,  to  describe  the  stress  intensity  (see  equation  (5-1)),  the  correction 
function  for  the  plate  with  hole  and  for  the  lug  as  Refs  61  and  72  works  out  as  follows: 


Y  (a  =  0)  = 


.  1.13  or  1.1. 


a  designates  the  stress  at  the  edge  of  the  hole.  The  validity  of  this  result  was 
max  °  J 

proved  theoretically  by  Neuber. 

Neuber  showed  in  Ref  76  for  the  crack  from  a  notch  through  a  boundary  value  analy¬ 
sis  a  ■*  0  that  the  factor  for  stress  concentration  results  from  the  multiplication 
of  the  factors  corresponding  to  the  macro  amd  micio  notch.  For  short  cracks 


\  =  Vr 

is  valid. 

Neuber  terms  the  factor  for  the  stress  concent  t  ..t  ion  of  the  macro  mo  toll  if)  and 

the  micro  notch  (=  crack)  a  a  is  equivalent  to  the  correction  function  Y  for  a 

K  K 

short  crack. 

Thus  a  precise  value  can  be  stated  for  the  correction  function  of  the  lug  at  the 
edge  of  the  hole. 

This  result  will  be  adopted  in  the  later  determination  of  the  stress  intensity 
for  the  various  lugs. 

6  DEF INITION  OF  PROBLEM 

As  already  mentioned  in  the  introduction,  lugs  are  essential  components  of  an 
airframe  for  which  proof  of  damage  tolerance  has  to  be  undertaken. 

Since  the  literature  does  not  contain  the  stress  intensity  solution  for  lugs  whii*’ 
is  required  for  proof  of  damage  tolerance,  the  problems  posed  in  the  following 
investigation  are: 

selection  of  a  suitable  method  of  determining  the  stress  intensity; 

determination  of  the  stress  intensity  as  a  function  of  the  crack  length  for 
various  forms  of  lug; 

showing  the  effect  of  essential  lug  parameters  on  sttess  intensity; 
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setting  up  a  complete  formula  for  calculation  of  the  stress  intensity  for  the 
lug,  allowing  for  essential  parameters; 

checking  and  evaluating  the  equation. 


Starting  point  in  laying  down  the  types  of  lug  to  be  investigated  are  the  diemsions  of 
lugs  in  use  in  aircraft  construction. 

If  the  lug  form  is  specified  in  a  simplified  manner  by  the  following  parameters 


radius  of  hole  R 
width  of  lug  W 
height  of  head  M  and 
thickness  t 

typical  lug  values  result  as  shown  in  Table  6.1. 

The  lug  dimensions  are  defined  in  Fig  6.1.  The  cheek  width  b  frequently  used 
later  is  calculated  from 


b 


W  -  2R 
2 


(6-1) 


L. 


The  geometry  of  the  lugs  to  be  examined  is  to  be  laid  down  in  such  a  manner  that  on  the 
one  hand  the  typical  lug  dimensions  in  Table  6.1  are  given  predominant  consideration 
and,  on  the  other  hand,  systematic  investigation  of  the  individual  parameters  in  regard 
to  their  effect  on  scress  intensity  is  possible. 

It  therefore  appears  logical  to  stipulate  a  lug  with  median  dimensions  and  median 
geometric  relations  by  means  of  the  values  quoted  in  Table  6.1.  For  the  lug  termed 
'standard'  below  the  following  mean  values  result 


Standard 


2R 

40  mm 

2R/W  = 

0.48 

2R/H  = 

0.9 

t 

1  5  mm. 

In  setting  down  the  lug  parameters  to  be  examined  one  cannot  rely  exclusively  on 
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the  few  investigations  quoted  in  the  literature  ’  .  In  these  tests  only  the  ratio 

between  hole  radius  and  lug  widtli  was  regarded  as  essential  parameter.  In  contrast, 

4-8 

the  numerous  fatigue  strength  tests  with  lugs  show  that  in  addition  to  parameter  2R/W 
the  ratio  of  2R/H  and  the  absolute  lug  size  are  also  significant  for  crack  life. 


In  a  comprehensive  lug  investigation  therefore  at  least  the  two  parameters  2R/W 
and  2R/H  should  be  considered.  Further  influencing  parameters  such  as  fit,  direction 
of  introduction  of  force,  etc,  are  not  dealt  with  in  this  paper.  This  is  reserved  for 
a  future  investigation  which  can  be  based  on  the  results  of  tiiis  paper. 

As  far  as  the  effect  of  size  is  concerned,  it  may  be  assumed  that  with  a  specific 
lug  geometry  at  the  same  nominal  stress  and  the  same  length  of  crack  the  speed  of  crack 
propagation  and  thus  also  the  stress  intensity  will  increase  as  the  size  of  the  lug 
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increases.  The  correction  function  which  very  probably  contains  the  dimensionless  crack 
length 

Y  =  f(a/b)  (with  b  as  cheek  width) 

as  free  parameter,  however,  ought  to  be  independent  of  the  lug  size.  This  is  shown  for 
instance  in  the  crack  model  "plate  strip  with  cracked  hole"  (see  equation  4-15)).  In 
order  to  demonstrate  this  on  the  lug  too  the  effect  of  size  will  also  be  taken  into 
account  in  the  investigation. 

If,  as  intended  in  this  investigation,  the  lug  is  regarded  as  a  plane  problem  the 
lug  thickness  t  is  not  important.  It  has  also  been  shown  in  fatigue  strength  tests 
that  within  a  technically  reasonable  range  the  thickness  of  the  lug  has  no  effect  on 
life. 

Table  6.2  lays  down  values  for  the  parameters  2R/W  ,  2R/H  and  the  effect  of 
size  expressed  by  2R  to  be  examined. 

7  INVESTIGATION  TO  DETERMINE  THE  STRESS  INTENSITY  FOR  Ll’GS  OF  VARYING  GEOMETRY 
7 .  1  Selection  of  a  suitable  method 

Section  3  presented  the  various  methods  which  can  be  used  to  determine  the  stress 
intensity  K  .  However,  very  few  are  suitable  for  application  to  the  lug.  Due  to  the 
complicated  boundary  conditions  the  analytical  methods  are  eliminated  1 rom  the  start. 

In  contrast,  the  method  of  finite  elements  has  proved  versatile  in  ’ts  use  for  a  variety 
of  crack  problems.  The  results  which  can  be  achieved  with  the  use  ot  this  numerical 
method  depend  decisively  on  how  well  the  structure  or  component  and  the  introduction  of 
force  can  be  idealised.  Idealising  the  lug  geometry  in  itself  presents  no  problems. 

But  when  the  introduction  of  force  and  its  effect  on  stress  distribution  and  stress 
intensity  is  regarded,  modelling  can  only  take  place  by  matching  to  a  known  stress  dis¬ 
tribution.  This  means  that  the  stress  distribution  for  a  lug  without  crack  would  tirst 
have  to  be  known.  For  this  case  the  selection  of  force  introduction  in  the  IT.  calcula¬ 
tion  would  have  to  be  such  as  to  achieve  the  best  possible  agreement  with  the  known 
stress  distribution.  Next  the  stress  intensity  lor  the  cracked  lug  can  be  determined. 

If  one  wanted  to  testify  on  the  accuracy  of  the  stress  intensity  determined  hv  the 
finite  element  method,  this  would  have  to  be  done  by  a  comparison  between  crack  propaga¬ 
tion  test  and  calculation.  For  these  reasons  the  finite  element  method  will  not  he  used 
for  the  systematic  examination  of  the  lug.  It  should,  however,  be  mentioned  hole  that 
further  work  with  the  finite  element  method  -  partial  results  are  already  available  'see 
Ref  82)  -  is  being  carried  out  especially  for  the  lug  and  other  typical  .tin  rail  compo¬ 
nents.  The  most  important  results  of  these  investigations  will  be  mentioned  insotar  as 
they  are  of  interest  here. 

Finally  there  is  nothing  left  but  experimental  methods  to  determine  the  stress 
intensity  for  the  lug. 


f 
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Of  the  various  experimental  methods  like 

crack  propagation 
comp  1 iance 
strain  gauges 
photoe las t i ci ty 

the  first  two  methods  are  particularly  of  interest.  While  the  strain  gauge  method  could 

lead  to  unacceptable  integration  errors  because  of  the  small  component  dimensions  and 

the  relatively  wide  strain  gauges,  photoelasticity  has  the  disadvantage  of  not  providing 

80 

particularly  accurate  absolute  values  .  A  further  drawback  of  the  last  two  methods  is 
that  the  stress  measurements  on  the  cracked  component  can  be  falsified  because  of 
yielding  in  the  area  of  the  crack  tip. 

Which  of  the  last  two  methods,  namely  the  compliance  method  described  in 
section  3.2.1  or  the  crack  propagation  method  described  in  section  3.2.2,  should  be  used 
for  the  lug  was  decided  on  the  basis  of  preliminary  tests  in  favour  of  the  crack  propaga¬ 
tion  method.  The  two  methods  were  compared  in  a  preliminary  programme  from  the  point  of 
view  of  scatter  of  test  results  and  cost  of  tests.  In  compliance  measurement  the 
difficulty  arose  that  after  each  crack  expansion  under  constant  amplitude  load  the  path 
measurement  equipment  had  to  be  readjusted  This  produced  deviations  with  the  compliance 
method  which  were  of  the  order  of  the  change  in  stiffness  to  be  expected.  This  was 
particularly  evident  in  the  area  of  short  cracks.  The  technical  problems  could  most 
probably  have  been  eliminated  by  appropriate  technical  effort.  But  this  would  appreciably 
increase  the  experimental  costs  which  were  almost  twice  as  high  for  the  compliance  method 
as  for  the  crack  propagation  method. 

The  essential  advantage  of  the  crack  propagation  method  can  be  recognised  particu¬ 
larly  from  an  engineering  point  of  view. 

Although  the  stress  intensity  can  be  determined  by  different  methods,  there  remains 
with  all  of  them,  insofar  as  they  are  approximation  methods,  the  uncertainty  in  regard 
to  the  accuracy  of  the  stress  intensity  solution.  Proof  can  only  be  provided  by 
additional  tests. 

This  proof  is  superfluous  if  the  stress  intensity  is  derived  from  crack  propagation 
test  results.  This  method  also  offers  the  advantage,  particularly  in  the  case  of  the 
lug,  that  the  introduction  of  force  does  not  have  to  be  idealised  as  in  many  other 
methods.  This  removes  a  not  inconsiderable  source  of  error. 

Another  reason  favouring  the  crack  propagation  method  is  that  a  considerable 
measure  of  experience  is  already  available  for  the  performance  and  evaluation  of  crack 
propagation  tests  from  which  this  investigation  can  benefit. 

Thus  the  technical  risk  can  be  reduced  and  the  quality  of  the  test  results 
improved . 

The  reliability  of  this  method  lias  also  been  demonstrated  in  an  investigation 
7  8 

carried  out  by  Hut h  to  determine  the  stress  intensity  in  bending  stressed  notches. 
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7 . 2  Specimen  and  test  set-up 

7.2.1  Dimensions  of  lugs 

In  the  fallowing  investigation  the  stress  intensity  is  to  be  determined  for  seven 
different  types  of  lug  by  crack  propagation  tests.  The  lugs  will  be  marked  1  to  7.  The 
lug  geometry  is  laid  down  in  Fig  7.1  by  the  hole  radius  R  ,  width  W  and  head 
height  H  . 

In  addition  the  thickness  is  designated  t  and  the  total  length  of  the  lug  L  . 

The  cheek  width  b  which  is  used  below  as  reference  quantity  is  given  by  b  =  (W  -  2R)/2  . 

The  dimensions  of  the  lugs  are  shown  in  Table  7.1. 

The  crack  length,  emanating  from  the  edge  of  the  hole,  is  termed  a  and  the  distance 
from  the  edge  of  the  hole  in  the  crack-free  component  x  . 

Figs  7.2  to  7. A  show  the  way  in  which  the  standard  lug  No. 2  is  varied.  Fig  7.2 
shows  the  different  lug  size  with  the  same  geometric  relations  (.represented  by  lugs  1,  2 
and  3)  . 

The  variation  of  cheek  width  with  the  same  hole  diameter  and  the  same  head  height 
is  shown  in  Fig  7.3,  with  lugs  4,  2  and  5.  The  radius  for  the  curvature  of  the  lug 
head  is  constant.  This  applies  also  to  the  variation  of  head  height  which  can  be  seen 
in  Fig  7.4  with  lugs  6,  2  and  7. 

7.2.2  Material 

Aluminium  alloy  7075  T7  351  from  Messrs,  t'egedur  was  used  for  the  lugs.  A  hOmm  thick 
plate  was  available  as  semi-finished  material.  The  lugs  were  taken  from  the  plate  in 
grain  direction. 

The  strength  and  chemical  composition  of  the  aluminium  alloy  is  shown  in  Tables  7.2 
and  7.3.  While  the  nominal  values  were  taken  from  the  material  specification  sheet,  the 
actual  values  come  from  workshop  tests  at  MMB,  CIVILE  No. 9337. 

For  tha  chemical  composition  the  actual  value  must  lie  within  the  nominal  value 
range.  For  strength  the  nominal  values  are  always  minimum  requirements  which  must  be 
reached  or  exceeded  by  the  actual  values. 

7.2.3  Test  set-up 

7.2. 3.1  Hydropulse  equipment 

A  hydropulse  equipment  from  Messrs.  Schenck  is  used  to  perform  the  crack  propaga¬ 
tion  tests. 

Fig  7.5  shows  a  schematic  presentation  of  the  entire  equipment  with  input, control 
and  recording  units. 

The  test  set-up  consists  of  a  base  plate,  two  vertical  columns  and  a  traverse. 

The  hydropulse  longitudinal  cylinder  is  incorporated  in  the  base  plate,  as  Fig  7.5 
shows.  The  specimen  clamp  is  placed  between  cylinder  and  traverse  and  is  mounted 
f lexib ly . 
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The  test  equipment  incorporates  an  inductive  displacement  transducer  on  the  cylinder 
and  a  dynamometer  beneath  the  ti averse  to  pick  up  the  test  values.  Additional  test  value 
receivers  can  be  fitted  to  the  test  object. 


The  test  signals  are  amplified  in  (?)  and  fed  to  the  recorder  and  the  controller  (? 
There  is  a  monitor  socket  (?)  between  test  value  pick-up  and  amplifier  (?) .  The 
controller  (?)  has  the  task  of  comparing  test  value  and  nominal  value  and  making  adjust¬ 
ments  through  the  servo  valves.  The  control  signal  is  amplified  in  (?)  . 

The  nominal  value  input  leads  via  the  monitor  socket  (?)  to  the  set  point  adjuster 
which  is  connected  to  the  generator  (?) .  The  generator  produces  periodic  sines, 
triangular  and  rectangular  functions.  The  static  nominal  value  is  preset  by  a  precision 
potentiometer  in  the  set  value  transmitter  (?)  and  controlled  by  the  digital  display  (?)  . 
The  set  value  is  passed  to  the  controller  (?)  . 

Wi  th  regard  to  ampl  i  tude  deviation  AA  the  manufacturer  quotes  the  following  value: 


AA  <  2%  up  to  f 


10  Hz 


AA  is  the  percentage  deviation  of  the  actual  value  from  the  set  value.  The  full 
amplitude  deviation  is  only  reached  at  a  limiting  frequency  of  f ^  '  30  Hz.  Since  the 
test  is  performed  at  a  considerably  lower  frequency  the  deviation  will  lie  far  below  22. 


7.2. 3.2  Installation  of  specimen  in  machine 

Fig  7.6  shows  the  installation  arrangement  for  the  lug.  The  end  of  the  lug 
opposite  to  the  hole  is  fixed  between  two  clamping  plates.  The  introduction  of  force 
into  the  hole  of  the  lug  takes  place  via  a  steel  pin  with  transition  fit  H7j6.  The 
steel  pin  is  lodged  on  both  sides  in  extremely  thick  steel  shackles  which  are  connected 
flexibly  with  the  traverse.  As  Fig  7.6  shows,  both  steel  shackles  have  observation  slits 
for  crack  propagation  measurement.  Since  this  entails  a  considerable  reduction  in  cross 
section  and  steel  shackles  had  to  be  made  thick  enough  for  the  hole  to  take  the  lug  pin 
not  to  penetrate  the  shackles  completely. 

Although  it  was  intended  to  examine  crack  propagation  in  one  direction  only,  obser¬ 
vation  slits  had  to  be  fitted  on  the  right  and  left  of  the  pin  and  in  both  shackles  for 
reasons  of  symmetry. 

7. 3  Test  procedure  and  test  results 

7.3.1  Test  procedure 

For  the  crack  propagation  test  a  crack  starter  notch  (see  Fig  7.7)  is  produced  by 
spark  erosion  at  the  edge  of  the  hole  of  the  lug. 

In  preliminary  tests  with  various  types  of  crack  starter  notch  the  so-called 
'V  notch'  proved  particularly  suitable  in  regard  to  the  generation  of  a  fatigue  crack 
with  the  straightest  linear  crack  front. 


A  similar  type  of  crack  starter  notch  but  with  an  acute  angle  is  already  used  as 
a  chevron  notch  in  the  compact  tensile  test.  The  shallow  angle  for  the  lug  tsee  Fig  7.7) 
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has  the  advantage  that  crack  propagation  measurement  can  be  started  at  relatively  short 
crack  lengths. 

Generation  of  a  fatigue  crack  with  approximately  straight  crack  front  offers 
following  advantages  for  determination  of  the  stress  intensity  from  crack  propagation 
measurements : 

The  crack  size  can  be  defined  by  a  single  parameter,  ie  the  crack  length  a 
(see  Fig  7.1),  in  contrast  to  the  semi-elliptical  crack. 

The  material  characteristic  da/dN  =  (AK^)  required  to  determine  the  stress 
intensity  is  established  with  standard  specimens  (central  crack  specimen, 
compact  tensile  specimen)  in  which  an  approximately  straight  crack  front  occurs. 
This  enables  logical  transfer  of  crack  propagation  characteristics. 

In  order  to  be  able  to  measure  crack  propagation  as  precisely  as  possible  the  lug  was 
provided  with  a  scale  along  the  crack  plane.  Using  a  microscope,  crack  growth  could  be 
read  off  to  an  accuracy  of  0. 1  mm. 

Crack  propagation  took  place  under  constant  amplitude  sinusoidal  loading  Isee 
Fig  7.8). 

? 

For  all  lugs  the  stress  was  =  98.1  N/mnT,  the  stress  a  being  related  to  the 
nett  cross  section. 


c  _ _ l _ 

N  (W  -  2R)t 


(7-1) 


where  P  indicates  the  force  transmitted  by  the  pin. 
The  stress  ratio  selected  was  R  =  0.  1  . 


R 


o 

u 


(7-2) 


The  test  frequency  was  10  Hz. 

Measurement  of  crack  propagation  started  at  an  initial  crack  length  of  a  =  2.5  mm. 

7.J.2  Tes  suits 

Three  crack  propagation  tests  were  carried  out  for  each  lug.  An  exception  was 
lug  No. I  for  which  an  additional  test  was  carried  out  because  of  somewhat  greater 
scatter.  The  test  results  for  the  seven  different  lugs  are  shown  in  Figs  7.9  to  7.15. 

The  number  of  load  cycles  to  fracture,  the  appropriate  mean  value  and  the  standard 
deviation  are  also  contained  in  the  relevant  figures.  The  small  scatter  for  crack 
propagation  tests  is  due  inter  alia  to  the  fact  that  the  specimens  were  taken  from  one 
batch  of  material.  A  mean  curve  was  plotted  through  the  scatter  band  of  the  test  values 
for  a  lug. 

The  mean  crack  propagation  curves  of  lugs  I,  2  and  !(=  effect  of  size)  are  shown 
in  Fig  7.16.  ft  will  be  seen  that  with  the  scale  enlargement  of  the  lug  the  critical 
crack  length  is  also  increased.  The  number  of  load  cycles  to  fracture,  on  the  other  hand, 

,.l 


is  not.  .if  t  tN'  t  ed .  This  iihmiis  that  tile  speed  of  ei.uk  propagation  also  rises  with  increas¬ 

ing  lug  sice.  Fig  '.17  shows  tin  mean  crack  propagation  curves  for  lugs  a,  and  r> 

(effect  of  cheek  width !  .  It  can  be  seen  that  with  increasing  cheek  width  the  critical 
crack  length  increases  but  the  number  of  load  cycles  Lo  fracture  declines.  The  mean 
crack  propagation  curves  lor  lugs  b,  7  and  .  can  be  seen  in  Fig  7.18.  The  increase  in 
head  height  produces  an  increase  ta  crack  propagation  life  while  the  critical  crack 
length  remains  the  same. 

7 .  -»  Discussion  of  test  results  with  reference  to  the  stress  distribution  in  the 

crack- free  lug 

To  determine  the  stress  intensity  for  an  infinite  plate  with  cracked  hole  and 
...  ,  ( )  ^ 

introduction  of  force  by  means  of  a  pin  Grandt  used  the  method  of  linear  superpos i t ion 
developed  by  Rice  (.see  section  4.7.1).  If  for  a  crack  model  with  a  sped  lie  introduction 
of  force  the  stress  intensity  and  the  crack  opening  displacement  as  a  function  of  the 
crack  length  V  =  flat  and  for  the  same  model  but  with  changed  introduction  of  force 
the  stress  distribution  arc  known  in  the  crack-fret?  model,  the  stress  intensity  for  the 
new  introduction  of  force  can  be  calculated  analytically  by  the  method  of  linear  super¬ 
position.  An  enquiry  into  this  method  in  regard  to  its  application  was  carried  out 
(see  Ref  79) . 

It  is  essential  for  this  method  that  the  stress  distribution  in  a  crack-free  state 
be  used  in  determining  the  crack  intensity.  Furthermore  Xeuber  proved  that  for  very 
small  cracks  the  stress  ot  the  macro  notch  enters  directly  into  the  correction  function  V 
(see  section  5).  It  must  therefore  be  assumed  that  the  stress  in  the  immediate  vicinity 
of  the  crack-tree  notch  has  a  vital  effect  on  the  stress  intensity  and  thus  on  crack 
propagation. 

In  Ref  80  stress  measurements  were  made  on  lugs  on  which  the  present  crack  propaga¬ 
tion  tests  were  later  carried  out.  1’hotoe  last  i  c  i  t  v  and  strain  gauges  were  used  as 
methods  to  determine  the  stress  distribution  in  the  fatigue  critical  lug  cross  section. 
The  object  of  this  stress  determination  was  to  develop  a  suitable  force  introduction 
model  tor  subsequent  examination  of  the  cracked  lug  with  the  aid  of  the  method  of  finite 
elements  (see  Kef  87).  The  stress  measurements  determined  (see  Ref  80)  lor  the  seven 
different  lugs  are  shown  in  Figs  7 .  1 l)  and  7 .  7  I  . 

If  these  results  are  compared  with  the  crack  propagation  measurements  the  following 
relationships  may  he  observed. 

Lugs  I,  7  and  1  in  Fig  7 . ] a  (  =  effect  of  size)  as  expected  show  the  same  stress 
iistribution  if  the  dimensionless  stress  o/z  with  -•  /I’tW  -  hit  is  plotted  over  the 
also  dimensionless  quantity  x/b  (=  distance  from  edge  ot  hole).  It  can  be  seen  from 
Fig  7.16  that  the  three  different  sized  lugs  exhibit  the  same  crack  propngnt ion  life. 

The  increase  in  critical  crack  length  connected  with  the  size  ot  lug  is  probable  due  to 
the  fact  that,  in  larger  lugs  the  boundary  el  feet  -  generally  leads  i  -  enhanced  increase 
in  stress  intensity  -  only  makes  itself  felt  at  longer  crack  lengths. 
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Fig  7.20  shows  the  stress  carve  for  lugs  4,  2  and  5.  It  is  noted  that  the  stress 
in  the  notch  area  increases  considerably  with  increa  :ng  cheek  width  b  .  As  can  be  seen 
from  Fig  7.17,  crack  propagation  life  declines  with  increasing  cheek  width  and  thus 
increasing  boundary  stress.  The  enlargement  of  the  critical  crack  length  connected  with 
the  increase  in  cheek  width  is  probably  also  attributable  to  the  boundary  effect  which 
only  makes  itself  felt  at  longer  crack  lengths. 

Fig  7.21  shows  that  the  stress  at  the  edge  of  the  hole  is  reduced  with  increasing 
head  height.  Comparison  with  Fig  7.18  shows  that  here  too  an  increase  in  stress  in  the 
immediate  notch  area  produces  a  drop  in  the  crack  propagation  life.  Since  the  cheek 
width  was  not  changed  in  this  case,  lugs  7,  2  and  6  yielded  the  same  critical  crack 
length. 

Comparison  between  crack  propagation  life  and  stress  in  the  crack-free  area  near 

the  notch  has  shown  that  the  crack  propagation  life  is  essentially  dependent  on  the 

stress.  Since  stresses  in  the  form  of  stress  concentration  factors  are  already  available 
75  81 

in  the  literature  ’  for  most  notches,  at  least  qualitative  statements  can  bo  made  with 
regard  to  the  crack  propagation  life  to  be  expected  if  the  shape  of  the  notch  is  changed. 

8  EVALUATION'  OF  THE  TEST  RESULTS  TO  DETERMINE  STRESS  INTENSITY 

8. 1  Premises  for  the  transfer  of  crack  propagation  characteristics  from  specimens 
to  structural  components 

Determination  of  the  stress  intensity  Kj  and  its  dimensionless  quantity  V  takes 
pLace  in  principle  by  the  method  described  in  section  3.2.2. 

This  method  is  based  on  the  fact  that  the  crack  propagation  rate  is  a  function  of 
the  stress  intensity  range 

da/dX  =  f(AK  )  . 

How  a  particular  value  of  is  composed,  ie  which  contribution  is  made  by  the  indi¬ 

vidual  terms  such  as  stress  o  ,  crack  length  a  or  correction  function  V  ,  is 
immaterial  provided  certain  conditions  are  followed: 

AKj  -  AnvwaY  . 

Transfer  of  crack  propagation  characteristic  determined  on  specimens  with  known  stress 
intensity  solution 

da /dN’  =  f(AK  ) 

to  components  is  permissible  if  in  determining  crack  propagation  for  the  specimen  and 
for  the  component  conformity  exists  in  the  following  points: 

Material 

Batch,  direction  of  specimen  extraction,  thickness  of  specimen. 

Environment 


Environment  medium,  temperature,  humidity. 
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Loading 

Frequency,  form  of  load  sequence,  stress  ratio  R  . 

Type  of  crack 

Formation  of  crack  front,  type  of  crack  opening. 

8 . 2  Determination  of  crack  propagation  characteristics 

To  determine  the  crack  propagation  characteristic  which  is  frequently  represented 
in  a  diagram  by  da/dX  over  2>F.^  ,  constant  amplitude  crack  propagation  tests  are 
performed  with  specimens  with  known  stress  intensity  solution.  As  a  rule  central  crack, 
compact  tensile  or  bending  spec i  — : ns  are  used  for  this  purpose.  A  central  crack  specimen 
with  dimensions  400  *  IbO  rnm  for  thickness  t  =  8  and  15  mm  is  used  here  for  the  crack 
propagation  tests.  For  a  material  thickness  of  t  =  20  mra  the  compact  tensile  specimen 
is  used.  The  compact  tensile  specimen  has  the  advantage  here  that  low  loads  are  required 
to  generate  crack  propagation  in  relatively  thick  specimens. 

Figs  8.1  and  8 . 1  show  the  two  specimens  with  relevant  stress  intensity  solution.  1 

In  order  to  ensure  the  transferability  of  the  crack  propagation  characteristic 
to  be  determined  to  lugs,  the  following  conditions  were  maintained  in  regard  to  the 
specimens  and  the  performance  of  the  test. 

Material 

AL  7U75T7  131  same  batch  as  for  lugs.  Specimen  extraction  in  L-T  direction 

l  =  longitudinal  side  of  specimen  in  direction  of  grain),  specimen  thickness 

t  =  8  and  15  mm  for  central  crack  specimen  and  t  =  20  mm  for  compact 

tensile  specimen. 

Environment 

Normal  laboratory  conditions. 

Loading 

Frequency  of  10  Hz,  sine-shaped  load  sequence,  stress  ratio  R  =  0. I  . 

-  Type  of  crack 

Crack  front  approximately  straight,  type  of  crack  opening  I  . 

The  test  set-up  described  in  section  7.2.3,  was  used  for  performance  of  tne  test 
but  with  suitable  specimen  clamping  device.  Crack  propagation  measurement  was  visual 
with  the  aid  of  a  microscope.  The  Lest  results  are  shown  in  Figs  8.3  and  8.4. 

The  test  results  show  that  scatter  of  the  crack  propagation  curves  -  for  the 
central  crack  specimens  three  were  carried  out  for  each  specimen  thickness  and  for  the 
compact  tensile  specimen  two  tests  only  -  is  relatively  small. 

A  direct  comparison  between  the  crack  propagation  measurements  on  compact  tensile 
specimens  and  central  crack  specimens  is  not  possible  since  the  two  types  differ  widely 
in  regard  to  stress  intensity.  Comparison  between  the  crack  propagation  measurements  on 
central  crack  specimens  shows,  however,  that  crack  propagation  life  for  the  8mm  specimen 
is  about  20%  longer  than  that  for  the  15mm  specimen.  Various  tests  have  shown  that 
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with  increasing  material  thickness  the  speed  of  crack  propagation  also  increases  and 
crack  propagation  life  declines. 


The  effect  of  material  thickness  on  the  speed  of  crack  propagation,  however,  is 

4  2 

shown  in  increased  measure  m  the  area  of  the  plane  stress  state  .  Since  the  material 
thicknesses  investigated  here  can  already  be  classified  in  the  transition  area  to  a 
plane  strain  state,  the  speed  of  crack  propagation  is  likely  to  increase  only  slightly 
with  increasing  specimen  thickness.  This  is  confirmed  if  the  crack  propagation  rate  is 
plotted  against  the  stress  intensity  range  (see  Fig  8.5). 

In  principle  the  crack  propagation  characteristic  can  be  obtained  by  determining 
the  slope 

m  =  .'.n/.'.N 


for  a  median  crack  propagation  curve  at  various  crack  lengths  and  the  appropriate  stress 
intensity  range  AK  .  The  quantities  required,  namely  the  range  of  stress  .V"  ,  crack 
length  a  and  correction  function  V  =  ffn/W)  are  known . 

Evaluation  of  the  crack  propagation  curves  becomes  simpler  and  more  accurate  if 
a  suitable,  constantly  differentiable  function  is  laid  through  the  measuring  points 
(=  regression  analysis).  The  Forman  equation  has  proved  a  most  suitable  regression 
function.  A  programme  was  developed  in  Kef  if  enabling  the  Forman  equation  to  be  matched 
to  crack  propagation  curves  determined  with  specimens  with  known  stress  intensity 
solution.  The  free  parameters  in  the  Forman  equation,  v'd  C  and  n^,  (see  equation  2.11) 
are  selected  so  that  the  least  squares  of  the  divergence  between  test  values  and  analytic 
function  show  a  minimum.  This  programme  lias  meanwhile  been  included  in  the  Handbook  for 
Structural  Calculation  (USB)  of  the  German  aircraft  industry  as  a  generally  recognised 
method . 


The  crack  propagation  characteristic  for  the  three  different  material  thicknesses 
determined  with  this  programme  shows  that  at  a  certain  stress  intensity  the  crack  pro¬ 
pagation  rate  increases  only  slightly  with  the  specimen  thickness.  Major  deviations  do 
not  occur  until  immediately  prior  to  reaching  the  critical  crack  length.  Fig  8.5  also 
includes  the  Forman  characteristics  for  the  relevant  material  thickness.  Here  it  should 

be  noted  that  the  Forman  characteristics  ,  n„  and  K.  differ  considerably  from 

F  F  c 

each  other  in  part.  If  specimens  of  different  thickness  exhibit  almost  the  same  crack 

propagation  behaviour  over  a  relatively  wide  range  one  would  expect  this  to  apply  to  the 

consistency  of  the  individual  characteristics.  However,  it  has  to  be  remembered  that 

C_  ,  n„  and  K  are  dependent  on  one  another.  But  if  the  K  value  depends  on  the 
F  1  c  c 

material  thickness,  the  C  and  n  values  must  change  to  the  effect  that  in  a  certain 

r  r 

AKj  range  the  Forman  equation  nevertheless  leads  to  similar  crack  propagation  behaviour. 


Determination  of  the  crack  propagation  rate  for  lugs  by  means  of  their  crack 


ion  curves 


Before  the  stress  intensity  or  its  dimensionless  quantity  V  can  be 

assigned  to  lugs  through  the  material  characteristic  da/dN  =  ffAK^)  (see  Fig  8.5)  it 


is  necessary  to  determine  the  crack  propagation  rate  da/dN  as  a  function  of  the  crack 
length  a  or  a/b  for  the  lug  crack  propagation  curves  (see  Figs  7.16  to  7.18). 

This  means  that  the  slope  of  the  crack  propagation  curves  must  be  determined  at 
different  crack  lengths. 

In  order  to  be  able  to  determine  the  slope  as  accurately  as  possible,  the  mean 
crack  propagation  curves  (see  Figs  7.16  to  7.18)  of  the  seven  lugs  are  approximated  by 
constant  differentiable  functions.  Crack  length  a  is  presented  as  a  function  of  the 
load  cycle  N  for  the  purpose. 

The  Forman  equation,  used  in  section  8.2  for  approximation  of  the  crack  propagation 
curves  of  the  central  crack  and  compact  tensile  specimen,  unfortunately  cannot  be  used 
here  since  the  stress  intensity  which  would  be  needed  is  not  yet  known  for  lugs  and 
still  has  to  be  determined.  An  attempt  is  therefore  made  to  approximate  the  lug  crack 
propagation  curves  by  polynomial  formulae.  In  doing  so  i t  has  been  found  that  the  crack 
propagation  curves  can  be  approximated  according  to  curve  piece  by  piece  by  two  or  three 
polynomials  of  the  second  order: 


■t 

a  =  <l]  +  q2N  +  q 


(8-1) 


q^_^  being  the  coefficients  to  be  determined.  Since  forevaluation  the  mean  lug  crack 
propagation  curves  are  already  used,  determination  of  the  coefficients  q.  was  carried 
out  not  by  'least  squares'  (regression  analysis)  but  by  introducing  three  pairs  of 
values  (a.,N.).  This  produces  three  equations  with  three  unknowns.  To  solve  the 
equations  pocket  calculator  HP-67  was  used  with  a  suitable  standard  programme.  This 
programme  handles  the  equation  system  in  the  form  of  a  vector  matrix  multiplication 


3  =  |Ml  x  q  .  (8-2) 

The  programme  first  forms  the  inverse  of  the  matrix  '  and  with  it  multiplies 

equation  (8-2) 


|M|  *  a  =  q 


(8-3) 


Equation  (8-3)  thus  solves  equation  (8-2)  by  the  unknown  coefficients  q.  . 

Coefficients  q  ^  ,  q.,  and  q^  for  the  individual  crack  propagation  ranges  are 
shown  in  Table  8.1  for  the  various  lug  crack  propagation  curves. 

The  crack  propagation  rate  is  obtained  by  di f ferent iat ing  equation  (8-1) 

da/dN  -  q.,  +  2q.(N  .  ( 8-., ) 

Since  the  polynomials  do  not  exhibit  equally  good  conformity  with  the  eraek  propagation 
curve  over  the  entire  sub-range  of  the  crack  propagation  curve,  da/dN  =  fia)  is 
determined  only  for  those  crack  lengths  or  load  cycles  for  which  the  polynomial 
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reproduces  the  crack  propagation  curve  as  accuratel'  as  possible.  i  tie-  crack  priipj^a!  i  mi 
rates  thus  established  are  shown  in  lips  H.n  to  8 .  8  tor  t  her  various  1  ups . 

Although  the  crack  propagation  rates  are  not  proportional  to  the  stress  intensity  - 
the  relationship  between  da/dN  and  Ak  .  >  ■  K  ^  appe.irs  linear  only  with  double  loga¬ 

rithmic  plotting  -  some  indications  can  already  he  gleaned  on  the  future  course  of  the 
stress  intensity.  All  curves  have  an  initially  steep  rise  shortly  after  the  start  of 
crack  propagation  in  common  which  in  all  tests  originated  at  an  initial  crack  length 
of  =  2.5  mm.  The  crack  propagation  ra  e  was  presumed  as  da/dN  =  0  at  the  initial 
crack  length  . 

Since  the  crack  length  is  relative  to  cheek  width  b  which  differs  for  the 
various  lugs  (see  Table  6.21,  the  points  at  da/dN  =  0  also  diverge  partially. 

The  slope  of  the  curves  declines  in  the  median  crack  length  range.  An  increased 
slope  is  then  noted  in  the  case  of  longer  cracks. 

As  Figs  8.6  to  8.8  show,  the  crack  propagation  rate  depends  to  a  very  large  extent 
on  all  the  influencing  factors  examined,  on  the  absolute  lug  size,  the  cheek  width 
2R/W  and  the  head  height  2R/1I  . 

The  effect  of  size,  seen  in  Fig  8.6,  causes  the  crack  propagation  rate  to  increase 
with  lug  size  at  the  same  crack  length  ratio  a/b  .  However  the  variation  in  cheek  width 
expressed  as  2R/W  has  a  considerable  effect  on  the  crack  propagation  rate.  As  shown 
in  Fig  8.7,  at  the  same  crack  length  ratio  a/b  the  crack  propagation  rate  increases 
considerably  with  cheek  width. 

On  the  other  hand,  the  effect  of  head  height  on  the  crack  propagation  rate  is 
less  pronounced.  Here  it  must  be  noted  that  with  increasing  head  height,  dei ined  by 
2R/H  ,  the  crack  propagation  rate  declines  for  a  certain  crack  length  ratio  t see  Fig  8.8). 

8.4  Determination  of  stress  intensity  and  its  non-dimensional  quantity  V  for 

the  various  lugs 

Using  the  material  characteristic  da/dN  =  l  ( ,',K  )  determined  in  section  8.2, 

appropriate  AK^  values  can  he  allocated  to  the  crack  propagation  rate  da/dN  determined 

on  lugs  for  different  crack  lengths  a/b  .  Basically  this  can  he  done  bv  extracting  I  torn 

Fig  8.5  the  appropriate  AK  ^  value  with  a  specific  da/dN  value  of  a  lug.  Since  for 

a  lug  the  stress  to  generate  the  crack  propagation  in  the  form  of  R  and  ’  and 

max 

the  crack  length  a  at  which  da/dN  was  determined  arc  known,  the  dimensionless  stress 
intensity  Y  ,  also  termed  correction  function,  can  he  caleiilalcd  as 


(1  -  R)  •  > ■ a. 

max  i 

This  procedure  can  he  carried  out  for  various  crack  lengths  lor  which  the  ct.uk  propaga¬ 
tion  rate  of  the  lug  is  known.  As  a  result  the  correction  function  lor  the  various  lugs 
can  be  presented  as  a  function  of  the  crack  length  a  or  its  dimensionless  quantity  a  b  . 
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The  use  of  the  material  characteristic  shown  in  Fig  8.5  for  assigning  the  stress 
intensity  to  the  crack  propagation  rate  of  the  lug,  however,  entails  the  risk  of  errors 
in  reading,  especially  since  the  plotting  is  double  logarithmic. 

This  source  of  error  can  be  removed  by  presenting  the  material  characteristic  in 
the  form  of  a  mathematical  function,  namely  the  Forman  equation.  It  has  already  been 
shown  in  section  8.2  how  the  Forman  equation  can  be  used  to  describe  crack  propagation 
curves  and  thus  also  the  crack  propagation  characteristic. 

The  connection  between  da/dN  and  AK^.  is  laid  down  in  the  Forman  equation  by  the 

characteristics  C  ,  n  and  K.  ,  quoted  in  Fig  8.5,  for  the  various  material 
r  r  C 

thicknesses : 

nF 

C  AK 

w*M»>lug  -  (r-RyK  -at:  •  <8-6 

c  1 


Since  the  equation  cannot  be  solved  according  to  Kj.  ,  must  be  determined 

iteratively. 


Here  again  a  standard  programme  of  the  I!l’fS7  is  used  to  determine  the  zero 
position 


0 


n 

CpAK  h 

R)K  -  AK, 
c  I 


da 

dN  ' 


(8-7) 


After  presetting  an  estimated  value  the  zero  position  is  calculated  by  the  following 
recurrence  formula 


X.  .  =  X.  -  f(X.) 

i  + 1  i  i 


(X. 

l 


v,) 


f (X . )  -  f(X.  .) 

1  1- 1 


(8-8) 


The  zero  position  or  the  corresponding  AK^.  value  is  determined  precisely  to  two  decimal 
places.  The  possible  error  is  less  than  1%. 

Expansion  of  this  programme  with  equation  (8-5)  produces  the  appropriate  correction 
function. 

The  course  of  correction  function  Y  depending  on  the  crack  length  a/b  for  the 
various  lugs  is  shown  in  Figs  8,9  to  8.11. 

It  is  revealed  that  the  Y  values  of  lugs  1,  2  and  3,  with  which  the  effect  of 
size  is  to  be  investigated,  lie  within  one  scatter  band  (see  Fig  8.9).  This  means  that 
despite  different  crack  propagation  behaviour  lugs  of  the  same  geometry  but  different 
sizes  exhibit  the  same  correction  function.  This  will  be  of  considerable  importance  in 
the  practical  application  of  the  correction  functions  determined  here. 

Based  on  the  considerations  in  section  5,  according  to  which  the  correction 
function  Y  for  a  -*■  0  followed  from  the  stress  concentration  factor  of  the  macro 
notch  and  a  factor  to  allow  for  the  load-free  surface,  the  experimentally  determined 
correction  function  was  interpolated  between  a  =  0  and  the  first  measuring  points. 
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distance  from  the  edge  of  the  hole.  The  stress  concentration  factor  and  thus  the  stress 
in  the  immediate  vicinity  of  the  edge  of  the  hole  increase  considerably  as  the  cheek 
width  b  gets  larger.  Enlargement  of  the  head  height,  in  contrast,  produces  a  reduction 
in  the  stress  concentration  factor. 

The  stress  concentration  factor  and  the  stress  in  the  immediate  vicinity  of  the 
edge  of  the  hole  therefore  have  a  great  effect  on  the  crack  propagation  life  and  govern 
the  size  of  the  correction  function  at  the  edge  of  the  hole,  which  was  pointed  out 
in  section  5. 

In  framing  a  formula  to  describe  the  correction  function  for  the  various  lugs 
the  stress  at  the  edge  of  the  hole  will  therefore  have  to  be  regarded  as  an  essential 
quantity. 

The  course  of  the  stress  at  a  greater  distance  from  the  edge  of  the  hole,  on  the 
other  hand,  does  not  allow  conclusions  to  be  drawn  on  crack  propagation  behaviour  or 
therefore  on  the  course  of  the  correction  function.  In  this  area  the  edge  effect  and 
the  rigidity  behaviour  of  the  cracked  component  will  be  significant  for  the  correction 
function.  The  formula  to  be  developed  will  therefore  have  to  rely  on  test  results  in 
this  area. 

9  DEVELOPMENT  CF  Adi  ANALYTICAL  FORMULA  TO  DESCRIBE  THE  CORRECTION  FUNCTION 

It  is  the  aim  of  the  following  reflections  not  to  describe  each  individual  experi¬ 
mentally  determined  correction  function  by  a  special  function,  such  as  for  instance  a 
polynomial,  but  to  cover  the  entire  test  results  in  a  complete  calculation  formula.  In 
principle  this  would  provide  a  means  of  describing  all  forms  of  lug  which  lie  within  or 
without  the  range  investigated.  This,  however,  presupposes  that  the  effect  ol  the 
individual  lug  parameters  on  the  correction  function  can  be  regarded  separately,  which 
has  been  confirmed  at  least  in  the  range  examined. 

Further  it  is  presupposed  that  the  connection  between  stress  intensity  and 
parameter  is  established  by  three  test  points  with  adequate  accuracy  at  least  within  the 
parameter  range  examined.  Outside  the  range  examined  any  statement  is  only  permissible 
to  a  limited  extent. 

Direct  evidence  of  these  assumptions  is  not  entirely  possible  since  an  infinite 

number  of  tests  would  be  required.  but  with  the  variation  of  lug  geometry  it.  may  lie 

assumed  with  adequate  safety  that  no  instability  occurs  in  regard  to  the  effect  on  the 

correction  function.  As  far  as  tin-  stress  at  the  edge  of  the  hole  is  concerned  -  the 

H  1 

stress  is  proportional  to  the  correction  function  there  -  Frocht  and  Hill  were  able 
to  show  that  the  effect  of  a  !  u>;  parameter  on  the  stress  at  the  edge  of  the  hole  exhibits 
a  steady  course. 

9 .  I  Computation  formula  for  _tje_  standard  lug 

aince  the  variation  of  tie  essentia]  I  parameters  emanated  from  lug  No..' 

(see  Table  7.1)  -  this  lug  was  therefore  termed  standard  lug,  -  a  formula  tot  the  •  or  re,  - 
tion  function  will  I  i  r  s l  he  developed  lor  tins. 
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As  shown  in  section  5,  the  following  applies  in  principle  to  cracks  from  notches 
with  load-free  surface: 


Y  =  K  1 .  12  for  a  -+  0  . 
t 


This  value  Y  at  the  edge  of  the  hole  is  used  for  all  lugs  as  so-called  'hanger'  of  the 
correction  function. 

Together  with  the  subsequent  decline  of  the  correction  function  -  similar  behaviour 
was  also  noted  in  the  stress  distribution  in  the  crack-free  lug  -  its  course  will  be 
determined  by  the  following  equation: 


1 . 1 2K  A 

v  _  _ t 

A  +  a/b 


(9-1) 


For  lugs  1,  2  and  3  according  to  Ref  81,  K  =  2.8  applies. 

Quantity  A  is  fixed  in  such  a  way  that  even  in  the  area  of  medium  and  longer 
crack  lengths  the  function  agrees  as  accurately  as  possible  with  the  correction  function 
determined  experimentally  in  Fig  8.9. 

If  equation  (9-1)  is  solved  in  respect  to  A 


1  .  1  2K  -  Y 

t 


(9-2) 


and  A  calculated  at  various  crack  lengths  a/b  with  the  appropriate  experimentally 
determined  Y  values,  it  is  established  that 

A  =  f (a/b) 


(see  Table  9.1  and  Fig  9.1). 

After  checking  various  functions  it  transpired  that  the  connection  between  A  and 
a/b  can  best  be  shown  by  an  exponential  function. 

The  free  parameters  in  this  exponential  function  are  calculated  with  the  help  of 
a  standard  programme  of  the  HP67  pocket  calculator  by  the  method  of  least  squares: 


A  =  0.026c'  .  (9 

A  comparison  between  the  experimentally  determined  values  for  the  correction 
function  and  the  values  calculated  to  equation  (9-1)  is  shown  in  Table  9.1.  The  mean 
error  over  the  entire  crack  length  range  is  approximately  1.8%. 

Since  there  is  no  size  effect  in  regard  to  the  correction  function  (see  Fig  8.9), 
equation  (9-1)  with  (9-3)  is  also  valid  for  lugs  I  and  3. 
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The  correction  function  can  thus  also  be  used  for  different  sized  lugs  with 
2R/W  =  0.48  and  2R/H  =  0.9. 

9 . 2  Extension  of  the  formula  to  take  into  account  the  effect  of  lug  cheek  width 

The  cheek  width  whose  effect  on  the  correction  function  was  examined  with  lugs 
No. 4  and  No. 5  (see  Fig  8.10)  is  characterised  by  the  ratio  2R/W  .  In  order  to  take  full 
account  of  this  effect  the  existing  formula  for  the  correction  function  (see 
equation  (9-1))  must  be  expanded  by  a  term  to  be  determined  empirically.  First  the  cor¬ 
rection  function  is  calculated  with  equation  (9-1)  for  lugs  No. 4  and  No. 5  at  various 
crack  lengths  a/b  .  The  following  values  are  used  (see  Ref  81)  for  the  stress  concentra¬ 
tion  factor  K 

t 

=  2.3  for  lug  No. 4 

K  =  3.85  for  lug  No. 5. 

The  values  are  shown  in  Table  9.2. 

Further  the  ratio  between  Y  /Y  is  determined  for  the  two  lugs  at 

test  calculation 

the  different  crack  lengths.  Y  means  the  correction  function  determined  in  the 

test 

test,  which  is  shown  in  Fig  8.10.  The  correction  function  calculated  to  equation  (9-1) 

is  termed  Y  ,  ,  .  .  The  individual  pairs  of  values  Y  and  a/b  and  the  ratio 

calculation  r 

Y 

test 

Y 

calculat ion 


are  also  shown  in  Table  9.2. 

Fig  9.2  shows  the  curve  of  Y  /Y  ,  ,  and  its  approximation  bv 

test  calculation  r 

equation  (9-4)  as  a  function  of  crack  length  a/b  for  lug  No. 4  (k^),  No. 5  (k^)  and 
for  No.  2  (k.,)  . 

For  lug  No. 2  the  additional  term  to  allow  for  2R/W  must  be  constant  over  the 
entire  crack  lengtii  range  and  assume  the  value  I. 

The  task  now  is  to  find  a  function  for  the  three  curves  k, ,  k„  and  k  (see 

4  2  5 

Fig  9.2)  which  fulfils  two  requirements.  First,  the  connection  between  k  and  a/b 
must  be  able  to  be  reproduced  with  adequate  accuracy.  Second,  the  function  may  only 
contain  one  free  parameter,  which  again  must  be  a  function  of  2R/W  .  The  following 
formula  presents  itself  as  a  suitable  function. 

r/a/b 

k  =  e  (4-4) 

r  being  a  function  of  2R/W  . 

In  order  to  be  able  to  describe  the  three  curves  in  Fig  4. 2  with  equation  (9-4) 
as  accurately  as  possible,  r  must  assume  the  following  values: 

r  =  0.25  for  lug  No. 4 

r  »  0  for  lug  No. 2 

r  -  -0.48  for  lug  No. 5. 
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If  k,  and  k.  are  calculated  with  these  r  values  and  compared  with  Y  / 

45  _  test 

Y  ,  ,  (see  Table  9.2)  the  result  is  for  lug  No. 4  a  mean  error  F  =  1.24 %  and 

calculation  _  ° 

for  lug  No. 5  a  mean  error  of  F  =  0.89%. 

The  relationship  between  r  and  2R/W  is  established  by  a  polynomial  second  degree. 
Since  the  three  r  values  are  reproduced  precisely  by  the  polynomial,  no  additional 
error  occurs  here. 

r  =  -  3.22  ♦  10.39(f)  -  7.67(f)~  .  (9-5) 

9 . 3  Extension  of  the  formula  to  take  account  of  the  effect  of  the  head  height  of  the  lug 

The  effect  of  the  head  height  of  the  lug,  character i sed  by  2R/H  ,  was  examined  by 
means  of  lugs  6  and  7. 

The  experimentally  determined  correction  functions  are  shown  in  Fig  8.11.  The 
values  for  Y  at  different  crack  lengths  are  shown  in  Table  9.3. 

If  the  correction  functions  are  calculated  to  equation  (9-1)  with  the  appropriate 
formulae  (see  Ref  81) 

K  =  2.6  for  lug  6 

K  =  3.0  for  lug  7 

and  the  ratio 

Y  /Y  ,  ,  . 

test/  calculation 


formed,  values  are  yielded  which,  irrespective  of  crack  length,  scatter  by  a  specific 
value  (see  Table  9.3). 

The  resultant  mean  value  is  used  as  additional  correction  U  to  allow  for  the 
effect  of  the  head  height. 

This  produces 


u  = 

0.97 

for 

lug 

6 

wi  th 

2R/H  =  0.7 

u  = 

1  .01 

for 

lug 

7 

wi  th 

2R/H  =  1.2 

u  = 

1.00 

for 

lug 

2 

with 

2R/H  =  0.9 

U  can  thus  be  shown  as  a  function  of  2R/H  . 

By  using  a  second  degree  polynomial  the  following  equation  is  yielded  for  U 


U 


0.72  + 


0.52 


(f-6) 


The  mean  deviation  or  the  mean  error  F  produced  by  using  U 

Y  over  the  entire  crack  length  range  is 

calculation 


instead  of 


'  test/ 


/I 


l 


,1 
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F  =  1.23 %  for  lug  6 
F  =  2.012  for  lug  7. 

If  equation  (9-b)  were  used  over  the  entire  crack  length  range  this  would  falsify  the 
correction  function  in  the  immediate  vicinity  of  the  edge  of  the  hole.  It  is  known 
that  at  the  edge  of  the  hole 


Y  =  1.12  for  a/b  =  0  . 

Therefore  1'  =  f(2R/H)  must  be  coupled  with  a  carrier  function  which  leads  to 

U  -*■  1  for  a/b  •*  0  . 


For  a/b  >  0  +  c  with  e  ^  1  the  carrier  function  must  not  have  any  effect  worth 
mentioning  on  U  . 

The  following  equation  suggests  itself  for  this  purpose: 


Q 


+  10 


-3 


10 


19-7) 


At  a/b  -+  0  function  Q  assumes  the  value  1  and  is  for  a/b  >  0  +  c  almost  identical 
with  U  .  The  error  is  less  than  0.06%  for  a/b  >  0.05  . 


As  may  be  gathered  from  the  data  in  Table  9.3  and  equation  (9-6),  the  effect  of 

the  lug  head  height  can  already  be  covered  to  a  very  large  extent  by  the  semi-empirical 

formula  (see  equation  (9-1)).  Minor  divergencies  between  the  correction  function  Y 

determined  experimentally  and  the  correction  function  Y  ,  ,  calculated  as 

calculation  . 

equation  (9-1)  Isee  Table  9.3)  are  allowed  for  by  the  function  Q  =  f  |l!(2R/H)  ^  (a/b)  J  . 

9 . 4  Stress  intensity  solution  for  the  lug 


if  the  functions  determined  in  sections  9.1  to  9.3  which  allow  for  the  effect  of 
the  different  lug  parameters  on  the  correction  function,  are  combined  they  produce  the 
following  solution  for  the  lug: 


K.  =  o  o  aY 

sum 


Y 

sum 
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A  +  a/b 
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II  £  +  10  * 
h 
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J.22  .  10.  »  (f ) 
0.72  •  0.52  (f)  - 
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This  solution  applies  to  the  unilaterally  cracked  lug.  The  loading  direction  must  point 
in  the  direction  of  the  principal  axis  of  the  lug.  The  crack  front  is  assumed  to  be 
approx i ina t ely  straight. 

Under  the  premises  laid  clown  in  section  '*  the  solution  is  not  only  applicable  with 
the  accuracy  demonstrated  in  sections  u.l  to  4.1  to  the  forms  of  lug  investigated  here. 

The  whole  purpose  ot  the  formula  is  to  cover  a  whole  spectrum  of  lug  forms.  It 
may  be  assumed  that  the  accuracy  of  the  stress  intensity  solution  for  lugs  within  the 
parameter  field  examined  should  not  alter  materially.  For  the  calculation  of  lugs  out¬ 
side  the  parameter  field  examined  a  .imilar  accuracy  can  be  assumed  if  the  'distance' 
is  not  too  great. 

The  types  of  lug  used  in  practice,  however,  should  be  covered  to  a  great  extent 
with  adequate  accuracy  by  the  formula  developed  here. 

1°  VERIFICATION  OF  THE  STRESS  INTENSITY  SOLUTION  FOR  THE  LUG  BY  MEANS  OK  EXPERIMENTALLY 

DETERMINED  CRACK  I’ROPACATIO.N  CURVES 

In  the  selection  of  a  suitable  method  of  determining  the  stress  intensity  for  a 
component,  section  7.1  pointed  to  the  advantage  of  the  crack  propagation  method,  namely 
the  possibility  of  checking  the  results  bv  means  of  crack  propagation  measurements. 

This  decisive  advantage  will  now  be  used  here  to  demonstrate  the  accuracy  of  the 
stress  intensity  solution  with  the  aid  of  the  c»-ack  propagation  measurements. 

Since  the  stress  intensity  solution  was  derived  directly  from  the  crack  propagation 
tests  there  should  theoretically  be  exact  agreement  between  the  crack  propagation 
investigation  using  the  stress  intensity  solution  (see  equation  (4-8)1  and  the  crack 
propagation  measurements. 

That  this  cannot  be  the  case  can  be  traced  in  the  main  to  two  causes  which  will  be 
dealt  with  in  more  detail  in  section  11  in  connection  with  considerations  of  error.  These 
causes  concern  the  unavoidable  material  scatter  and  the  divergencies  in  the  approximat i on 
of  the  experimentally  determined  correction  function  by  analytical  functions. 

The  crack  propagation  calculation  was  performed  with  the  Forman  equation 
(equation  (11-11))  by  numerical  integration,  using  a  standard  programme  ot  the  Hl’f>7  pocket 
calculator. 

This  programme  is  used  for  calculation  of  the  integral  ol  any  constantly  differen¬ 
tiable  functions  within  given  interval  limits. 

The  result  is  equal  to  the  area  which  the  function  encloses  within  the  limits 


with  the  abscissa. 
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Starting  from  the  initial  crack  length  the  number  of  load  cycles  N  to  crack 

length  a  is  obtained  by  integration  of  the  reciprocal  value  of  the  Forman  equation. 

The  variable  in  the  Forman  equation  is  the  stress  intensity  factor  range  AK^  which 
with  a  given  lug  geometry  and  external  load  only  remains  a  function  of  the  crack 
length  a  . 

To  perform  the  calculation  the  standard  programme  required  for  integration  is 
first  read  into  the  calculator  store  by  means  of  a  magnetic  card.  The  input  of  the 
Forman  equation  follows  in  which  in  AK^  the  entire  lug  correction  function  is 
contained.  After  input  of  the  number  of  integration  steps,  the  initial  crack  length 
and  the  crack  length  up  to  which  the  number  of  load  cycles  is  to  be  computed,  the  result 
yielded  is  the  appropriate  load  cycle  N 

(2 

The  necessary  input  data  are  shown  in  Table  10.1  for  the  seven  lugs. 

As  in  the  test,  the  stress  for  the  crack  propagation  calculation  for  the  seven 
lugs  was  cJq  =  98.1  N/rtutT,  and  the  stress  ratio  R  =  0.1. 

Figs  10.1  to  10.3  show  a  comparison  between  the  mean  experimentally  determined 
crack  propagation  curves  and  the  load  cycles  calculated  to  Forman  for  several  crack 
lengths.  Relatively  close  agreement  between  test  and  calculation  can  be  noted.  The 
deviation  between  calculation  and  test  is  approximately  10"  on  average.  The  greatest 
deviations  arise  for  lug  3  with  1 7%  and  lug  7  with  22".  Experience  has  shown  that 
deviations  of  the  order  found  here  are  within  the  range  of  material  scatter. 

1 1  INVEST IGAT ION  INTO  POSSIBLE  SOURCES  OF  ERROR 

As  can  be  seen  from  Figs  7.9  to  7.15,  the  crack  propagation  curves  determined  for 
lugs  are  subject  to  a  certain  scatter.  Since  the  correction  function  for  the  lugs  was 
derived  from  these  crack  propagation  curves  this  too  suffers  from  a  corresponding  error. 
Although  the  mean  crack  propagation  curves  were  used  for  the  evaluation,  the  correction 
function  cannot  be  regarded  as  the  mean  value  of  the  basic  totality  due  to  the  fact 
that  it  is  not  a  case  here  of  a  basic  totality  but  a  random  sample  of  limited  extent. 

A  further  error  of  approximately  the  same  siz.e  can  be  traced  to  the  scatter  behaviour 
of  the  material  characteristic  required  for  evaluation. 

As  far  as  the  correction  function  in  the  form  of  a  semi-empirical  formula  is 
concerned  (see  equation  <‘*-8),  this  too  is  subject  to  a  certain  error. 
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First,  possibilities  of  error  which  can  he  traced  to  the  scatter  of  experimentally 
determined  crack  propagation  curves  will  be  examined  more  closely. 

11.1  Scatter  of  crack  propagation  curves  and  its  ef f ect  on  the  correction  function 

Before  the  scatter  occurring  in  crack  propagation  tests  is  determined  it  must  first 
be  defined.  One  can,  for  instance,  understand  by  scatter  the  scatter  of  crack  length 
at  a  specific  load  cycle  and  in  reverse  the  scatter  ..1  the  load  cycle  at  a  specific  crack 
length.  But  one  can  also  regard  the  scatter  of  the  speed  of  crack  propagation  as  a  random 
event . 

No  matter  which  criterion  is  taken  as  random  event,  the  scatter  connected  with  crack 
propagation  measurements  can  be  traced  to  various  causes.  In  the  main  these  are  the 
material,  the  test  specimen,  the  test  set-up  and  the  reliability  of  the  reading  of  the 
person  performing  the  test.  It  is  to  be  assumed  that  these  sources  of  error  are  not 
constant  and  do  not  occur  at  the  same  intensity  during  the  whole  of  the  experiment. 

It  therefore  seems  to  make  sense  to  examine  only  the  effect  of  the  various  sources 
of  error  on  the  crack  propagation  life,  ;Y  the  sum  of  their  scatter  behaviour. 

This  simplifies  consideration  of  error  materially  and  enables  it  to  be  redu  ed 
to  the  following  question. 

If  crack  propagation  life,  :'c  the  load  cycle  to  fracture,  is  subject  to  certain 
scatter  behaviour,  what  is  the  error  for  the  correction  function  which  may  be  expected? 

For  the  rest,  the  consideration  of  error  will  be  undertaken  bv  the  basic  rules  of 

.  .  81  .  .  ... 
statistics  on  the  assumption  of  a  Gaussian  normal  distribution. 

Since  ihe  life  of  the  individual  lugs  is  of  different  lengths  but  it  may  be  assumed 
that  the  scatter  behaviour  is  independent  of  lug  geometry  to  the  greatest  possible  extent, 
the  variance  coefficient  will  be  used  to  describe  scatter  behaviour. 


(11-1) 


I,  signifies  the  mean  value  of  fatigue  life  and  s  its  standard  deviation. 
The  following  variance  coefficients  are  yielded  for  the  seven  lugs. 


0.103  0.027  0.0h3  0.04b  0.018  0.083  0.028 


Further  investigation  will  allow  for  the  mean  as  well  as  the  maximum  variance  coefficient, 


V  0.1 0  )  . 

max 
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Furthermore,  investigation  of  error  will  not  be  carried  out  on  every  lug  but  only  on  lug  2 
as  representative  (=  standard  lug). 

To  determine  the  relevant  variance  coefficient  for  the  correction  function  the 
crack  propagation  life  will  first  be  calculated  with  the  aid  of  the  Forman  equation 
(equation  (8-6))  with  the  mean  value  of  V  as  equation  (9-8).  All  crack  propagation 
calculations  are  performed  by  numerical  integration  of  the  reciprocal  Forman  equation 
with  the  HP67  pocket  calculator. 

The  crack  propagation  calculation  is  repeated,  assuming  different  standard  devia¬ 
tions  for  the  correction  function  until  the  variance  coefficient  of  the  calculated  life 
values  agrees  with  those  observed  in  the  test. 

In  regard  to  the  variance  coefficient  this  yields  .e  following  connection  between 
life  and  correction  function. 

Variation  coefficient 


Life 

Correction  function 

Maximum  value 

0.103 

0.039 

Mean  value 

0.05 

0.014 

Shown  in  a  probability  diagram  the  following  relationships  are  produced  between  the 
scatter  behaviour  of  the  crack  propagation  curve  and  the  correction  function  (see 
Fig  11.11). 

If  it  is  assumed  that  the  mean  velue  of  the  random  sample  in  relation  to  the  basic 
totality  does  not  occur  at  50%  (=  100%  life),  but  at  10%  or  90%  value  the  error 
arising  in  the  correction  function  is  2%  or  5%  respectively  for  V 

max 

In  experimental  evidence  of  life  it  is  customary  if  few  test  values  are  available 

to  allot  to  them  a  failure  probability  of  90%.  The  actual  mean  value  is  then  smaller 

.  .  84'  . 

according  to  the  standard  deviation  .  Since  at  least  three  tests  per  lug  were  performed 
in  this  investigation  allocation  of  the  mean  value  at  the  90%  value  is  fairly  conservative. 
For  this  reason  the  error  will  only  be  derived  from  the  mean  variance  coefficient.  The 
error  of  the  experimentally  determined  correction  function  i;  therefore  assumed  at  2"’. 

If  it  is  further  assumed  that  the  material  characteristic,  which  was  also  determined 
by  crack  propagation  tests,  may  be  subject  to  the  same  error  then  the  possible  error  of 
the  correction  function  is  increased. 

Thus  the  total  error  of  the  experimentally  determined  correction  function  may  reach 
2.83%.  According  to  Fig  1!.]  this  would  result  in  an  error  of  approximately  12%  for 
fatigue  life. 


The  divergence  between  the  experimentally  determined  correction  functions  (see 

Figs  8.9  to  8.11)  and  the  semi -empirical  formula  as  equation  (9-8)  can  be  shown  by  simple 

28  29 

error  observation.  It  is  customary  in  fracture  mechanics  ’  to  state  the  accuracy  of 
a  regression  equation  over  the  entire  crack  length  range  by  the  fol lowing  process: 
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Y  -  Y  .  ' 

test.  calculation 


H  1-2) 


F  is  the  mean  percentage  deviation.  n  is  the  number  of  points  at  which  a  comparison  is 
carried  out.  Based  on  this  consideration,  the  following  mean  deviations  are  produced  for 
the  correction  functions  of  the  seven  lugs  (see  Tables  9 . 1  to  9.3). 

l  ug _ 1,2,  ) _ _ 5 _ b _ 7__  . 

V"  1-8  l./'i  li.su  1..'  <  2.0 

The  divergence  or  accuracy  of  equation  (10-8)  over  the  entire  crack  length  range  is 
therefore  between  0.89?  and  2. O'7. 

If  the  crack  propagation  is  calculated  after  Forman  taking  account  of  the  correction 
function,  the  correction  function  can  be  subject  to  an  error  of  approximately  i.57  - 
errors  from  material  scatter  and  regression  -  at  maximum. 

The  individual  errors  were  added  according  to  the  following  formula 


■  w 


(11-3) 


with  F. 

l 

and  F 


individual  errors 
the  sum  of  all  errors. 


An  error  of  approximately  1 57  can  thereby  arise  for  crack  propagation  life.  It  is 
known  from  experience  that  errors  oi  this  order  are  in  the  area  of  natural  material 
scatter  in  crack  propagation  tests. 

If  crack  propagation  tests  are  performed  on  specimens  which  either  stem  from 

different  material  batches  or  from  semi-finished  materials  from  different  manufacturers 

8  5 

the  error  can  amount  to  up  to  100? 

Re-calculation  of  the  experimentally  determined  crack  propagation  curves  for  the 
seven  lugs  (see  Figs  10.1  to  10.3)  has  shown  that  the  divergence  between  test  and  calcula¬ 
tion  amounted  to  10%  on  average. 

12  CONCLUDING  OBSERVATIONS 

Higher  requirements  in  safety  and  increased  exploitation  of  the  capacity  of  material 
to  bear  stress  have  led  to  the  extension  of  fracture  mechanics  to  the  most  varied 
technical  areas. 

An  essential  constituent  of  the  calculation  principles  for  prediction  of  crack 
propagation  and  the  critical  crack  length  is  the  stress  intensity  solution. 

It  was  the  object  of  this  paper  to  determine  for  an  important  structural  component , 
namely  the  lug,  stress  intensity  as  a  function  of  the  typical  lug  parameters  and  present 
it  in  the  form  of  a  complete  analytical  formula. 

First,  the  stress  intensity  was  described  starting  from  linear  elastic  principles 
of  strength  and  its  importance  for  determination  of  crack  propagation  and  critical  crack 
length  demonstrated. 
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A  review  ol  the  best  known  methods  ol  determining  the  stress  intensity  followed 
together  with  a  diseussion  of  their  app l ieab t 1 i tv  to  the  problem  of  the  lug. 

This  revealed  that  determination  o!  the  stress  intensity  bv  means  of  crack  propaga¬ 
tion  curves  is  a  particularly  suitable  method  t or  the  lug. 

At  the  time  this  work  was  carried  out  there  was  no  comprehensive  investigation  in 
existence  on  stress  intensity  for  lugs,  especially  on  the  effect  of  the  various  lug 
parameters. 

Starting  with  an  infinite  plate  with  hole,  through  the  strip  with  loaded  and 
unloaded  hole,  tor  which  known  stress  intensity  solutions  were  available,  the  lug  was 
therefore  approached  step  by  step. 

In  a  subsequent  experimental  investigation  of  the  lug  iL  was  possible  to  demonstrate 
the  connections  between  the  essential  lug  parameters,  crack  propagation  behaviour  and 
stress  distribution  in  the  crack-free  lug,  in  addition  to  determining  the  stress  intensity 

The  set  aim  of  the  task  was  finally  achieved  by  succeeding  in  presenting  the  stress 
intensity  or  its  dimensionless  quantity  Y  as  a  function  of  the  essential  lug  parameters, 
such  as  cheek  width,  head  height  and  absolute  lug  size  in  a  complete  semi-empirical 
formula.  The  accuracy  of  the  formula  was  demonstrated  by  re-calculating  the  crack  pro¬ 
pagation  curves  for  the  various  lugs.  The  difference  between  test  and  calculation  was 
22%  maximum  and  approximately  10%  on  average. 

Strictly  speaking,  the  present  stress  intensity  solution  is  valid  only  for  the 
unilaterally  cracked  lug  with  an  approximately  straight  crack  front. 

In  practice  cracks  of  this  type  occur  predominantly  when  caused  by  fretting 
corrosion.  However,  cracks  are  frequently  found  in  lugs  which  emanate  from  the  edge 
of  the  hole.  The  crack  front  then  assumes  the  form  of  a  quarter  ellipse  or  quarter  circle 

There  is  as  yet  no  stress  intensity  solution  specially  for  this  case.  If  one  were 
to  use  the  stress  intensity  solution  for  a  crack  with  straight  crack  front  here  the  crack 
propagation  life  thus  calculated  would  lie  considerably  on  the  safe  side. 

The  following  procedure  should  lead  to  a  materially  more  accurate  result. 

If  the  strip  with  open  hole  is  considered,  there  is  a  stress  intensity  solution  for 
a  crack  with  straight  track  front  and  for  a  semi-elliptical  edge  crack1  . 

It  can  be  assumed  in  lirst  approximation  that  the  relationship  of  the  individual 
stress  intensity  solutions  is  transferable  to  the  lug. 

Nevertheless,  it  remains  for  future  investigations  on  the  lug  lo  expand  the  stress 
intensity  solution  for  the  crack  with  straight  or  approximately  straight  crack  Iron!  in 
order  to  be  able  to  take  account  of  the  corner  crack  and  the  effect  of  different 
directions  of  Load  application. 
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EFFECT  OF  LUC  PARAMETERS  ON  STRESS  DISTRIBUTION.  CRACK  PROPAGATION  AND  ON  THE 


CORRECTION  FUNCTION 


Parameter 


Effect  of 


Effect  of 
cheek  width 


Stress  distribution 


Same  stress  curve  (see 
Fig  7.19) 


Considerable  increase 
in  stress  at  edge  ot 
hole  with  increasing, 
cheek  width;  reversal 
with  increasing  dis¬ 
tance  from  edge  of 
hole  a/b  N  0.  I  S 
(see  Fig  7.70; 


Crack  propagation 


Same  crack  prop,  life; 
increase  in  crack  pro¬ 
pagation  rate  with  lug 
size  (see  Fig  7.16) 


Same  correction  func¬ 
tion  (see  Fig  8.9) 


Reduction  in  crack  prop.  Ris 
life  with  increasing  fun 
cheek  width;  increase  in  ing 
crack  growth  rate  with  j r in 
!  increasing  cheek  width  i  Kev 
!  -  see  F  i  g  7 .  I  7 )  1  a/i' 


i  with  mere.) 
'k  wi  1  f  h  i  n 


r  liijir  a/b  •  ().  ->. 

Reversal  in  range 
a/i'  ••  1 1 . 6  i  see 
Fi g  8.10; 


Effect  o  f 
head  height 


Rise  in  stress  at  edge 
of  hole  with  decreas¬ 
ing  head  height; 
rove rsa  I  at  a/h  ■  0 . 0  . 
(see  Fig  7  . I  ) 


(eduction  in  cr.nl;  prop,  j  Ri  s*.  in  c-  -rrec  t  i  < 
tile  will:  .le.Teas  i  nr  !  function  with 

h>-  1  * i  height;  higher  decreasing  head 

•  i  n  K  p  i  op  .  ra  i  ■  a  I  !  he  i  gill  i  nara  1  lei 

low.  •-  head  heigh'  j  .1  i  sp  1  acem.  a  t  i  -n 

he.-  :  ig  ' .  1  8  )  '  Fig  8.  '.I .) 


DATA  TO  DETERMINE  111!,  I.MFlKloAl  'M  ANTI  TV  A 


Lug  No. 6  Lug  No. 5  Lug  No. 


Table  9.2 


DATA  TO  DETERMINE  THE  EFFECT  OF  WIDTH  2R/W 


Y  test 

Y  calculation 

Y  test 


Y  calculation 
k4 

Error  Z 


Y  calculation 
k5 

Error  % 


0.1 

0.2 

1  .92 

1.64 

1  .  74 

1.44 

1  . 1 

1.14 

1  .08 

1.12 

1  .82 

1.75 

2.5 

1.95 

2.92 

2.41 

0.86 

0.81 

0.86 

0.81 

0 

0 

0.5 

0.6 

O 

CO 

O 

3.41  2.92  2.41  2.18  2.07  2. 031  2.04  I  2.08  1  2.141  2.22  I  2.31 


0.71  0.68  0.66  0.65  0.65  i  0.65 


0.7)  I  0.69  I  0.67  I  0.65  I  0.63  I  0.62 
0 


Table  9.3 

DATA  TO  DETERMINE  THE  EFFECT  OF  LUC  HEAD  HEIGHT 


0.05  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8 


2.16  1 .88  I  I .62 

2.31  1 .97  I  1.631  1.47 


v  .  ,  ,0.94  0.95  1  0.99  I  0.98 

Y  calculation 

^6 

Error  Z 


Y  test 


3.19  12.11  I  2. 


2.7  2.4  2.0 


Y  calculation 

according  2.66  2.27  1.88 

G 1  (9.1) 


Y  test 


Y  calculation 


Error  X 


1.01  1.05  1.06 

1.01  - 

0  3.8 


0.96  0.96  0.97  0.97  0.97  1  0.99 


0  0 


1 .58  I  1.57  1 .62  1.6 
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INDEX  OF  MOST  IMPORTANT  ABBREVIATIONS 


General  geometric 

quanti ties 

a 

ram 

length  of  crack 

ao 

ram 

initial  length  of  crack 

a 

c 

ram 

critical  length  of  crack 

U,  V 

ram 

displacement  in  x  and  v  direction 

W 

mm 

width  of  specimen 

t 

mm 

thickness  of  specimen 

R 

mm 

radius  of  hole 

b 

ram 

cheek  width  of  lug 

H 

mm 

head  height  of  lug 

L 

ram 

length  of  lug 

V  \ 

stress  concentration  factor 

y 

correction  function 

Quantity  related 

to  material 

E 

> 

.N/mm~ 

modulus  of  elasticity 

V 

Poisson’ s  rat io 

V 

'  e 

N/mm 

specific  surface  energy 

Y 

N  /mm 

specific  plastic  deformation  action 

p 

c 

N/mm 

stiffness 

c_1 

mm/N 

compliance 

S’  P  ;  CF’ 

nF 

crack  propagation  characteristics  for  Paris  or  Forman  equation 

S 

N/mm 

crack  resistance 

Klc 

N  /  mm-* !  - 

critical  stress  intensity  in  plane  strain  state 

K 

c 

N/mm-*/- 

critical  stress  intensity  in  plane  stress  state 

Quanti  t  ies 

related  to  load 

P 

N 

force 

a 

N/mm“ 

gross  stress 

aN 

N/mm*' 

nominal  stress 

ao 

N  /ram^ 

2 

maximum  stress  of  a  stress  cycle 

a 

N/mm 

minimum  stress  of  a  stress  cycle 

a 

N/mm^ 

stress  amplitude 

S 

AKt 

KfO 

N/mmi/2 

stress  intensity  for  crack  opening  mode  I 

N/mmi/2 

stress  intensity  range 

N/mm3/2 

stress  intensity  for  an  infinite  plate  with  central  crack 

G 

N/mm 

crack  extension  force 

I 
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Fig  2.1  Stress  sequence  with 
constant  amplitude 


Fig  2.2  Crack  propagation  rate  as  a  function 
of  the  stress  intensity  ranae 
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Fig  3.2  Idealisation  of  compact  tensile 
specimen  by  finite  elements 


Key: 

Vergrosserung  =  magnification 
Rissspitze  =  crack  tip 


Fig  3.1  Idealisation  of  a  plate  with  central 
crack  by  finite  elements 


Fig  3.3  Isoparametric  element 


Figs  3. 4-3. 7 


Fig  3. 


Fig  3.4  Simulation  of  stress  singularity 
by  an  isoparametric  element 


5  Idealisation  of  the  crack  tip 
with  the  aid  of  isoparametric 
elements 


Fig  3.6  Idealisation  of  the  crack  tip 
with  the  aid  of  conventional 
elements 


Key: 

Tuglast  =  tensile  load 

Verlangerung  =  extension 


Fig  3.7  Experimental  determination  of 
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Figs  3.8-3.10 
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Fig  3.9  Plot  of  isochromatics 
in  the  area  of  a 
sharp  notch 


Fig  3.8  Procedure  in  determining  stress 
intensity  with  the  aid  of 
the  crack  propagation  method 
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Fig  3.10  Determination  of  stress  intensity 
for  an  edge  crack  specimen 


Figs  4. 1-4. 4 


Fig  4. 
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Fig  4.1  Loaded  and  unloaded  hole  with  crack  in 
infinite  and  finite  plate 


2  Uni-  and  bilateral  crack  from  unloaded 
hole  in  an  infinite  plate 


Fig  4.3  Open  hole  with  crack  under 
load  p(x) 


Key: 

Unbelastete  Bohrung 
Spiel 

Dimensionslose  Risslange 


=  unloaded  hole 
=  clearance 
=  dimensionless  crack 


length 


Fig  4.4  Curve  of  stress  intensity  for  an  open  hole 
and  for  a  loaded  hole  with  different  pin  fit 


Figs  4.9-4.11 


Key: 

(un)belastete  Bohrung 
mit  zwei  Rissen 
Dimensionslose  Risslange 


=  (un) loaded  hole 
=  with  two  cracks 
=  dimensionless  crack  length 


Fig  4.9  Comparison  of  correction  functions  for  the  loaded 
and  unloaded  hole 


Fig  4.10  Stress  distribution  in  the 
non-era eked  lug 


Fig  4.11  Correction  function  for  a  lug 
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Fig  4.12  Effect  of  cheek  width  R/W  on  the  plot  of  the 
function  for  the  uni-  and  bilaterally  cracked 
the  method  of  linear  superposition 
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Fig  4.13  Finite  element  model  for  the 


a)  Punktlatt 


b)  Cotmui-Lattvaftailung 


c)  Laitvartailung  nach  Cotinusquadrat 


Key: 

(a)  concentrated  load  at  centre  of  hole 

(b)  cosine  load  distribution 

(c)  cosine  square  load  distribution 


Fig  4.14  Different  means  of  idealisation  of 
introduction  of  force  into  the  lug 
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Key: 

R'echnung  (untere  etc) 
Lastspiele 


=  calculation  (lower  and  upper 
scatter  limit) 

=  stress  cycles 


Fig  4.15  Comparison  between  test  and  calculation 
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Crack  model 


Sketch 


Description 


Infinite  plate  with 
O'  crack  from  open  hole 


Infinite  plate  with 
crack  from  hole  and 
force  introduction 
by  means  of  a  pin 


Plate  strip  with 
open  hole 


Plate  strip  with 
hole  and  force 
introduction  by 
a  pin 


Me  t  hod 


Analytical  stress 
function  conformal 
mapping 


Analytical  linear 
superposition  with 
aid  of  'weight 
function' 


Numerical  solution 
of  stress  function 
series  (boundary 
collocation) 


Numerical  solution 
of  stress  function 


Superposition,  super¬ 
imposing  known 
solutions 


Analytical,  linear 
superposition  with 
aid  of  'weight 
function ' 


Superposition,  super¬ 
imposing  known 
solut ions 


De  scribed 


section : 


Figs  5.1&6.1 


Re  fe  rerue 

No. 


4.1  Bowie  (6 1 ) 


4.2.2  Grandt  (64) 


4.3.1  Newman  (63) 


4.3.2. 1  Newman  (63) 


4. 3. 2. 2  Broek  (69) 


Impel  1 i zzeri 
and  Rich 
(72) 


Liu  and 
Kan  (73) 


Fig  5.1  Survey  of  different  crack  models  and  methods  to  determine  stress 
stress  intensity 
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Fig  6.1  Definition  of  notations 
for  the  lug 
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Figs  7. 1-7. 4 


Fig  7.1 


Fig  7 


Laying  down  lug  geometry  by 
typical  parameters 


.2  Variation  of  lug  by  scale 
enlargement 


Fig  7.3  Variation  of  lug  width  W 


Fig  7.4  Variation  of  head  height  H 


it  or 


Ragistnerung 


Sane  31  Grundausrustung*  Moduln 


Reqelstracke 


0«  .u>gf aph 


Terminal 


$offw«rtgeb<ar 


i — ; — - 1 

|L_, 

|l  ! 

HyJ'opun 

P'OfCiirechner 

®i 


a 


Sedienteia 


Key: 

Reg’ strierung 

X-Y-Schreiber 

Rechner-Terminal 

Sol lwertgeber 

Functionsgerierator 

Prozessrechner 

Bedienfeld 

Geschlossener  Regelkreis 
Regel strecke 

Priifaufbau  mit  Hydropul s- 
Langszyl inder 
DruckOl versorgung 
Hydraul ikaggregat 


=  recording 
=  X-Y  recorder 

-  computer  terminal 

=  setting  transmitter 
=  function  generator 
=  process  computer 
=  operating  field 
=  closed  control  circuit 

-  control  section 

-  test  set-up  with  Hydropulse 

longitudinal  cylinder 
=  pressure  oil  supply 
=  hydraulic  control  unit 


Fig  7.5  Plan  of  hydropulse  equipment  -  Series  31:  Basic  equipment  module 


f- 


Figs  7. 68.7.7 


Key: 

Einspannlasche  =  shackle 

Bolzen  =  pin 

Sehschl itz  =  inspection  slit 

Riss  =  crack 

Fig  7.6  Mounting  device  for  the  lug 
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in 


Fig  7.7  Crack  starter  notch  at  the  edge  of  the  hole 
in  the  lug 


Key: 

"Spannung  =  stress 
Zeit  =  time 

Fig  7.8  Constant  amplitude  load  sequence 
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Key: 

Lastspiel zahl  bei  Bruch  =  load  cycle  on  fracture 
Mittelwert  -  mr>an  value 

Standardabweichung  =  standard  deviation 


Fig  7.9  Crack  propagation  investigation  of  lug  1 
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Figs  7.10&7.11 
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Fig  7.10  Crack  propagation  investigation  of  lug  2 


Fig  7.11  Crack  propagation  investigation  of  lug  3 


Figs  7.12&7.13 
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Fig  7.12  Crack  propagation  investigation  of  lug  4 
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Fig  7.13  Crack  propagation  investigation  of  lug  5 
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Figs  7.1487.15 


Fig  7.15  Crack  propagation  investigation  of  lug  7 
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Figs  7.18&7.19 


Fig  7.18  Mean  crack  propagation  curves  for  lugs  2,  6  and  7 


Key: 

Bezogener  Abstand  vom  =  relative  distance  from 
Bohrungsrand  edge  of  hole 


Fig  7.19  Stress  distribution  for  lugs  1,  2  and  3 
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Key: 

Dicke  =  thickness 

Gul tigkei tsbereich  =  area  of  validity 

Genauigkeit  =  accuracy 


Fig  8.1  Central  crack  specimen  for  determination 
of  crack  propagation  characteri stic 
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Fig  8.2  Compact  tensile  specimen  for  determination 
of  crack  propagation  characteristic 


Fig  8.3  Crack  propagation  measurements  on  compact 
tensile  specimens 
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Figs  8.4&8,B 
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Fig  8.4  Crack  propagation  measurements  on  central 
crack  specimens 
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Fig  8.5  Crack  propagation  characteristic  for  aluminium  alloy  7075T7351 
with  differing  specimen  thickness 
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Fig  9. 
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Fig  9.1  Relation  between  the  empirically 
determined  quantity  A  and  crack 
length 


2  Relationship  between  correction  function  determined 
experimental ly  and  to  equation  (9-1)  as  a  function 
of  crack  length  for  lugs  4  and  5,  as  well  as 
approximation  with  equation  (9-4) 


Figs  10.1-10.3 
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Fig  10.1 
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Fig  10.2 
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Fig  10.3 

Key: 

Versuchsergebris  =  test  result 
Gerechnet  =  calculated 


Figs  10.1-3  Comparison  between  crack  propagation  curves  determined 
experimentally  and  calculated  as  Forman,  using 
equation  (9-8) 
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